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Abstract 

In this paper, we investigate thermal effects of the Jaynes-Cummings model 
(JCM) at finite temperature with a perturbative approach. We assume a single 
two-level atom and a single cavity mode to be initially in the thermal equilibrium 
state and the thermal coherent state, respectively, at a certain finite low temper- 
ature. Describing this system with Thermo Field Dynamics formalism, we obtain 
a low-temperature expansion of the atomic population inversion in a systematic 
manner. Letting the system evolve in time with the JCM Hamiltonian, we examine 
thermal effects of the collapse and the revival of the Rabi oscillations by means 
of the third-order perturbation theory under the low-temperature limit, that is to 
say, using the low-temperature expansion up to the third order terms. From an 
intuitive discussion, we can expect that the period of the revival of the Rabi oscil- 
lations becomes longer as the temperature rises. Numerical results obtained with 
the perturbation theory reproduce well this temperature dependence of the period. 
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1 Introduction 



The Jaynes-Cummings model (JCM) was originally proposed for describing the sponta- 
neous emission in a semi-classical manner by Jaynes and Cummings during the 1960s 
[U [21 El E]- The JCM consists of a single two- level atom and a single cavity mode of 
the electromagnetic field. The JCM interaction between the atom and the cavity mode 
is obtained by the rotating wave approximation, so that each photon creation causes an 
atomic de-excitation and each photon annihilation causes an atomic excitation. The JCM 
is a soluble quantum mechanical model. Moreover, it is simple enough for expressing the 
basic and most important characteristics of the matter-radiation interaction. Because of 
these distinct advantages, both theoretical and experimental researchers in the field of 
quantum optics have been studying the JCM eagerly for decades. 

If we initially prepare the atom in the ground state and the cavity mode in the coherent 
state, the JCM shows the periodic spontaneous collapse and the revival of the Rabi 
oscillations during its time-evolution [HI El [7J El [9] . This phenomenon was experimentally 
demonstrated in the 1980s [10J. We can regard this phenomenon as a direct evidence for 
discreteness of energy states of photons. Thus, the JCM has a fully quantum property, 
which cannot be explained by semi-classical physics. 

In this paper, we investigate the non-dissipative JCM at finite low temperature. We 
assume that the atom and the cavity mode are initially in the thermal equilibrium state 
and the thermal coherent state, respectively, at a certain finite low temperature j3[= 
l/{k^T)]. Thus, we can describe an initial probability distribution of quantum states 
of the system with the canonical ensemble. Moreover, we assume the time-evolution 
of the system is governed by a unitary operator generated with the JCM Hamiltonian. 
This implies that the system does not suffer from dissipation and its time-evolution is 
reversible. 

The Thermo Field Dynamics (TFD) formalism was developed by Takahashi and 
Umezawa for dealing with phenomena in isolated systems during the 1970s [TT j [T2" l IT5 1 IT3] . 
The TFD formalism has a wide range of applications in equilibrium situations of closed 
systems, as well. In this paper, we think about applications of the TFD to non-dissipative 
closed systems. The TFD formalism helps us to be more successful than conventional for- 
malisms do as follows: calculating an expectation value of a pure state created by the 
TFD mechanics, we can obtain a thermal average in statical mechanics. In return for this 
benefit, the TFD formalism requires us to let the original Hilbert space double in size for 
the tensor product. Then, the TFD formalism induces thermal-like noises in two mode 
squeezed states generated by the thermal Bogoliubov transformation. 

Describing the JCM at finite low temperature with the TFD formalism, we can write 
down a thermal average (an expectation value) of an observable as a series expansion 
containing powers of 8 ((3) in a systematic manner. [The explicit form of the function 6(f3) 
is given by 9((3) = arctanh[exp(— f3hu/2)}, where a; is a frequency of the cavity mode. We 
note that 9(f3) — > +0 as f3 — > +oo.] We call this series the low-temperature expansion. 
This prescription is a new key point of this paper as compared with the other past works. 

Strictly speaking, because we introduce a finite temperature into both the atom and 
the cavity field, the perturbation theory has to include two parameters of the temperature. 
In this paper, we let 0(/3) and 6((3) denote the parameters of the temperature for the 
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atom and the cavity field, respectively. [The explicit form of the function Q(/3) is given 
by 0(/3) = arctan[exp(— /3fiu /2)], where u is a transition frequency of the two-level 
atom. We note that ©(/?) — > +0 as — > +00.] However, the Hilbert space of the atom 
is two-dimensional, so that we can solve the problem concerning the atom exactly. On 
the other hand, the dimension of the Hilbert space of the cavity field is infinite, so that 
we cannot give a rigorous treatment for the problem concerning the cavity field, which 
interacts with the atom. This is the reason why the low-temperature expansion contains 
powers of 9(0), but not 0(/3). 

As mentioned before, in this paper, we initially prepare the cavity mode in the thermal 
coherent state, which is proposed by Barnett, Knight, Mann and Revzen as a natural 
extension of the zero-temperature ordinary coherent state according to the TFD [T5l [16] . 
Then, we investigate the time-evolution of the atomic population inversion using the 
low-temperature expansion up to the third order terms in 9(/3). We regard this low- 
temperature expansion as the third order perturbation theory. 

To take another approach, we give an intuitive discussion and obtain the following 
expectation: the period of the revival of the Rabi oscillations becomes longer as the 
temperature rises, where the temperature is finite but low enough and varies in the neigh- 
borhood of the absolute zero. Numerical results obtained by the third order perturbation 
theory reproduce well this temperature dependence of the period. 

Here, we write about related works. The thermal JCM without dissipation is inves- 
tigated in Refs. [TTj, [HI HH (201 [211 122] • In these references, the cavity field is prepared 
initially in a thermal equilibrium state, whose statistical behavior is given by the Bose- 
Einstein distribution. Then, the system consisting of the atom and the cavity mode is 
assumed to evolve with a unitary operator generated by the JCM Hamiltonian, so that 
its time-evolution is reversible. However, in these works, thermal coherent states are not 
considered to be the initial states of the cavity field. 

The thermal JCM with dissipation is investigated in Refs. [231 122 [251 1251 123 • In 
Refs. [231 [2H [23 [2S] , the cavity damping in the JCM is discussed. In these works, the 
equation of motion for the density operator p includes the term nL ir (p), which describes 
the irreversible motion caused by the cavity damping. A typical form of L- n (p) is given by 
(2apa) — pa) a — a)ap) in Refs. [2U [2S], where a) and a are the creation and annihilation 
operators of the cavity photons, respectively. On the other hand, in Ref. [21], the phase 
damping in the JCM is discussed. In this work, the equation of motion for the density 
operator p includes the term (— j[H, [H, p]]), which describes the phase damping. To 
obtain this term, we assume the system to interact with a heat-bath environment, namely 
an infinite set of harmonic oscillators. Considering the Liouville-von Neumann equation 
for the system and the heat bath, and applying the Markovian approximation to it, we 
can derive the term of the phase damping. 

Some researchers attempt to extend the JCM Hamiltonian according to the TFD 
formalism [15*1 [221 [2H]- However, in these works, they introduce the temperature only 
into the cavity field and let the atom maintain the temperature at absolute zero. In this 
paper, we explain that such extended JCM Hamiltonians cannot be genuine ones based 
on the TFD formalism. 

The organization of this paper is as follows: In Sec. [21 we give a brief review of 
the original JCM. In Sec. [31 we give a brief review of the TFD. In Sec. HI we give an 
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intuitive discussion about the thermal effects found in the period of the revival of the Rabi 
oscillations. In Sec. [5j we formulate the perturbation theory under the low-temperature 
limit. In Sec. El we compare our TFD formalism applied to the JCM with the Liouville- 
von Neumann equation. Although the equilibrium TFD is equivalent to the Liouville- 
von Neumann equation in principle, we find that the TFD formalism is more effective 
and easier than the Liouville-von Neumann equation for deriving the low-temperature 
expansion. In Secs.[7J[8]and[9l we calculate the zero-th, first and second order perturbation 
corrections, respectively. In Sec. [TUJ we give numerical results of calculations. In Sec. dU 
we consider thermal effects of the counter-rotating terms. In Sec. [121 we give a brief 
discussion. In Appendix |AJ we give details of calculations of the third order perturbation 
correction. 

The evaluation of the perturbation corrections described in Sees. El El [9] and Ap- 
pendix [AJ we make use of techniques for calculations developed in Ref. [22] ■ Thus, we 
can regard this paper as a sequel of Ref. [22J . 

2 A brief review of the JCM 

In this section, we review the JCM briefly. The Hamiltonian of the original JCM is 
expressed in the form, 

H = —tu a z + hua^a + hn,(a + a + cr_a^), (1) 
2 

where h = h/2ir, a± = {l/2)(o x ± ia y ) and [a,a^] = 1. The Pauli matrices (crj for 
i — x,y, z) are operators acting on quantum states of the single atom in the cavity. The 
creation and annihilation operators a) and a act on quantum states of the single cavity 
mode. The transition frequency of the two-level atom is given by uq. The frequency of 
the single cavity mode is given by u. From now on, for simplicity, we assume k to be real. 
We can divide the Hamiltonian H given by Eq. (TjQ) into two parts as follows: 

H = %{C X + C 2 ), 
d = cj{^<y z + a 1 a), 

C 2 = K,(a + a + cr_a f ) - ~^-^z, (2) 

where Acu = u — ojq. We can confirm [Ci,C2] = 0. Moreover, we can diagonalize C\ at 
ease. 

Because of the above facts, we take the following interaction picture for expressing 
the time-evolution of the quantum state: First, we write down the state vector of the 
atom and the cavity photons in the Schrodinger picture as \ips(t)). Second, assuming 
IV'i(O)) = iV's(O)}) we define the state vector in the interaction picture as \ipi(t)) = 
exp(iCit)\if)$(t)) . Thus, we can describe the time-evolution as \ipi{t)) = U(t)\ipi(0)), 
where U(t) = exp(—iC2t). 

We give eigenstates of the two-level atom and the cavity mode in the following forms: 
First, we describe the ground state and the excited state of the atom as two-component 
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vectors, 

|0>a=|O>a=( J), |e>A = |l>A=(j), (3) 

where |<?)a(= |0)a) an d I 6 )a(= |1)a) are eigenvectors of a z , and their corresponding 
eigenvalues are (—1) and 1, respectively. The index A stands for the atom. We can 
regard \i)a for % e {0, 1} as the number states of the fermions. Second, we write down the 
number states of the cavity photons as |n)p = (1/Vn\)(af) n \0)p for n = 0,1,2,.... The 
index P stands for the photons. 

We describe the unitary operator for the time-evolution U (t) as the 2x2 matrix, 

U(t)=e Xp [-it(-* u / 2 rA = { Ul1 UW ), (4) 
y no) Au/2 J' y U i MOO J 

where 
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, r~, . . . Au ^sin(V at a + c + l\n\t) 
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and 



u o = cos (y at a + c\n\t) — i— — -\fc , (5) 

| Aw | VaTa + c 

«-<£)". (0) 

Because we take the basis vectors {|1)a, |0)a}, the indices i,j G {1, 0} for Uij are arranged 
in descending order. That is to say, we take the index '1' for representing the first row 
and the first column of the 2x2 matrix U(t), and we take the index '0' for representing 
the second row and the second column of the 2x2 matrix U(t). 

After these preparations, the probability for detecting the ground state of the atom 
|(?)a at the time t is given by 

P 9 (t) = \\A(gMt))\\ 2 - (7) 

Moreover, the atomic population inversion is given by 

<(7,(0> = (^(0k.W*)> = i-2P ff (0- (8) 

As a particular case, we consider the initial state to be in |^i(0)) = |<7)a|c*)p, where 
\a)p represents the coherent state, 

\a)p = exp(o;a^ — a*a)|0)p 

oo n 

= e-'^E-^Wr, (9) 
«=o v n\ 
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and a is an arbitrary complex number. From now on, for simplicity, we always let the 
parameter a characterizing the coherent state |a)p be real. Then, we obtain P g (t) in the 
form, 



P 9 (t) 



(0 1) 


I Mil 




( i4p ) 




\ M 01 


MOO / 





'(a|MooMoo|a)i 



e _ Q 2 [ cos 2 (x/ri + cklt) + c 

rr'n n! 



2/i 



sm"(\/n + c|/c|i). 



(10) 



n=0 



n + c 

and ( imp shows the collapse and the 



The atomic population inversion given by Eqs. 
revival of the Rabi oscillations. 

Here, we examine the time scale of the initial collapse and the period of the revival of 
the Rabi oscillations [29]. Learning from experience, we know this phenomenon becomes 
more distinct as a 2 increases. Thus, we assume a 2 3> 1. We rewrite the index of the 
summation n as n = a 2 + 8n in Eq. ( ITOi) . Then, because of the Poisson distribution, the 
major contribution for the summation in Eq. ( jTOj) comes from \Sn\ < a 2 . Moreover, for 
simplicity, we assume c <C a 2 , and we neglect the term [c/(n + c)] sin 2 (>/ra + c|ac|£) in the 
right-hand side of Eq. (fit)]) . 

From the above discussions and Eqs. flSj) and ([TP]) , writing y/n + c ~ (a 2 + n)/(2|a|), 
we obtain the following approximation: 



M*)> 



00 2n 1 

^ n 1 ^i(a 2 +n)|K|t/|a] + e ~i(a 2 +n)| K |t/|a|l 



E 

n=0 



n! 2' 



- exp[a 2 (cos — 1)] cos(|a||/c|i + a 2 sin -LL 
a a 



)■ 



In Eq. (ITT]). exp[a 2 (cos(|fi;|t/|Q!|) — 1)] represents the amplitude envelope of the wave, and 
cos[|a||«|t + a 2 sin(|/c|i/|a|)] represents the Rabi oscillations. Therefore, we can estimate 
the time scale of the initial collapse and the period of the revival of the Rabi oscillations at 
and 27t|o;|/|k| around, respectively. Moreover, paying attention to a 2 sin(|/«|t/|a|) ~ 
I a I \n\t for a 2 ^> 1, we can estimate the period of the Rabi oscillations at about 7r/(|o!||«;|). 



3 A brief review of the TFD 

In this section, we give a brief review of the TFD developed by Takahashi and Umezawa 
[TT1 [I2l IT3l [14] . The TFD is a method for describing the quantum mechanics at finite 
temperature. Using this formalism, we can describe the statistical average of an observ- 
able at finite temperature as a pure state expectation value. Thus, if we take the TFD 
formalism, we do not need to deal with a mixed state, which is a statistical ensemble of 
pure states at finite temperature. 

In return for the above advantage, the TFD lets us introduce the so-called tilde par- 
ticles corresponding to the ordinary particles. Then, we understand that the ordinary 
particles and the tilde particles represent the dynamical degree of freedom and the ther- 
mal degree of freedom, respectively. Thus, to construct the TFD formalism, we introduce 
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a fictitious Hilbert space T-L corresponding to an original Hilbert space H and handle 
quantum mechanics on % ® H. 

In this section, according to the TFD formalism, we define the thermal vacua of bosons 
and fermions. Moreover, we discuss the thermal coherent state proposed by Barnett, 
Knight, Mann and Revzen [T5l [T6] . After these preparations, we rewrite the Hamiltonian 
of the JCM according to the TFD. 

First, we consider the TFD formalism for describing the system of the bosons. We 
define the ordinary Hilbert space 

H b :{|0)b,|1)b,|2) b ,...}, (12) 

and its corresponding tilde space, 

H b :{|0)b,|1)b,|2) b ,...}. (13) 

Then, the TFD formalism for the bosons is defined on the following space: 

{\n) B ® \fh) B E H B ® U B ■ n,m G {0,1,2,..}}. (14) 

We write the creation and annihilation operators on the Hilbert space l-i B as a) and a, 
respectively. Moreover, we write the creation and annihilation operators on the Hilbert 
space "Hb as a) and a, respectively. Then, we assume the commutation relations, 

[a,aJ] = [a, a)] = 1, [a, a] = [a,a^] = 0. (15) 

Next, we introduce the temperature (3 = l/(k B T) as follows: 

U B (9) = exp[z0(/3)G B ], (16) 

Gb = 2(aa — aW), (17) 

cosh^(/3) = [1 -exp(-/3e)]- 1/2 , 

sinh#(/3) = [exp(/3e) - 1]- 1/2 , (18) 

where e = hu. We note that the relations G B = G B , U B (9) = U B ^O) and U B (-9) = U B (6) 
hold. In Eqs. f|T6|) and (fT7|) . to emphasize that Gb and U B {9) are operators acting on both 
1-L B and H B , we put an accent (a hat) on them. Moreover, we pay attention to the fact 
that 9((3)(> 0) is real and 9((3) — > +0 as (3 — > +oo in Eq. ( ITS]) . The index B appearing in 
Gb and U B {9) stands for the boson. 

Because of introducing the temperature, the creation and annihilation operators de- 
fined on "Hb and "Hb are transformed as follows: 

a -> a{9) = U B (9)aU B (9) = cosh9(/3)a - smh9(f3)a\ 

a a{9) = U B (9)aU ] B (9) = cosh9(/3)a - smh9((3)al (19) 

The transformation given by Eq. f lT9|) is called the Bogoliubov transformation for two 
mode squeezed states. From Eqs. (|15p and f)TTJ|) . we can derive the commutation relations, 

[a(9), tf(6)] = [a(9), rf{9)] = 1, [a{9), a(9)] = [a{9), a\9)} = 0. (20) 
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Thus, we can regard a(0) and a(0) as quasi-particles at the temperature (3. 
Here, we define the zero-temperature vacuum as 

|0,0) B = |0)b®|0) b gH b 8)^b. (21) 

Furthermore, referring to Sec. 2.2.3 of Ref. [H], we define the thermal vacuum in the 
form, 

|O(0)) B = C>b(0)|O,6>b 

= exp(-lncosh^)exp[(tanh^)a t a t ]|0,0) B . (22) 

From Eqs. ( !T9|) and ( )22|) . we obtain 

a(0)|O(0)>B = a(0)|O(0)> B = O. (23) 

Thus, we can consider the thermal vacuum |O(0))b to be a vacuum for the quasi-particles, 
which are represented by a(0) and a(0). 

Second, we consider the TFD formalism for describing the system of the fermions. We 
define the ordinary Hilbert space, 

«p:{|0) F ,|l) P }, (24) 

and its corresponding tilde space, 

H F :{|0) F ,|1) F }. (25) 

Then, the TFD formalism for the fermions is defined on the following space: 

{|n) F ® |m) F eH F ®H F : n,m G {0,1}}. (26) 

We write the creation and annihilation operators on the Hilbert space H F as and c, 
respectively. Moreover, we write the creation and annihilation operators on the Hilbert 
space Hp as $ and c, respectively. Then, we assume the anti-commutation relations, 

{c, ct} = {c, c f } = 1, {c, c} = {c, ct} = 0. (27) 

Next, we introduce the temperature /3 as follows: 

U F (9) =exp[i0(/9)G F ], (28) 

Gp = i(cc — c^c^), (29) 

cos#(/3) = [l + exp(-/3e)]- 1/2 , 

sin 0(0) = exp(-/3e/2)[l + exp(-/3e)]~ 1/2 . (30) 

We note that the relations G F = G F , f> F (0) = U ¥ \0) and U F (-0) = U F (6) hold. More- 
over, we pay attention to the fact that 0(/3)(> 0) is real and 0(/3) — > +0 as (3 — > +oo in 
Eq. ( 130|) . f/ F (0) in Eq. ( 128|) and G F in Eq. ( )29|) are again of the same forms as Ub(0) in 
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Eq. f|T6|) and Gb in Eq. (fTTj) . respectively. The index F appearing in U F (9) and G F stands 
for the fermion. 

Because of introducing the temperature, the creation and annihilation operators de- 
fined on Hp and T-L F are transformed as follows: 

c c(6) = U F (9)cU F (9) = cos 9(/3)c + sin 9(f3)c\ 

c -> c(0) = U F (6)cUl(9) = cos9((3)c-sm9((3)cl (31) 

From Eqs. ( 127)) and ( 13 ip . we can derive the anti-commutation relations, 

{c(9),J(9)} = {5(9), d\9)} = 1, {c(9),c(9)} = W)^{e)} = 0. (32) 

Thus, we can regard c{9) and c{9) as quasi-particles at the temperature j3. 
Here, we define the zero-temperature vacuum as 

|o,o) F = |o) F ® |6) F e n F ®n F . (33) 

Furthermore, referring to Sec. 2.4.2 of Ref. [H], we define the thermal vacuum in the 
form, 

|O(0)) F = U F {9)\0,0) F 

= [cos 9 + (sin^)c t c t ]|0,0) F . (34) 

From Eqs. ( 13T)) and ( 134)) . we obtain 

c(9)\0(9)) F = ~c(9)\0(9)) F = 0. (35) 

Thus, we can consider the thermal vacuum \0(9)) F to be a vacuum for the quasi-particles, 
which are represented by c{9) and c{9). 

Third, we consider the thermal coherent state. After the above preparations, Barnett, 
Knight, Mann and Revzen define the thermal coherent state as follows 



\a, f, 6) B = exp[aa\9) + fa\9) - a*a{9) - ^a{9)\ |O(0)) B , (36) 

where a{9) and a{9) are given by Eq. ( ITU)) . 

However, the TFD formalism requires all state vectors to be invariant under the tilde 
conjugation, which is given by 

(XY)~ = XY, 
(Xt)~ = x\ 

[X)~ = aX, (37) 

where 

a={\ J° S ° n ) . , (38) 
1—1 (rermion) 

X and Y are arbitrary operators defined on 7i B and 7-L F , and and £ 2 are arbitrary 
complex numbers. For example, the thermal vacua are obviously invariant under the tilde 
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conjugation, that is to say, |O(0)) B ~ = |O(0)) B , b(O(0)|~ = b(O(0)|, \0(6)) f ~ = \0(9)) F 
and F (O(0)|~ = F (0(6»)|. Moreover, we can show (G B )~ = -G B , (U B (0))~ = Ub(0), 
(G F )~ = -G F and (U F (9))~ = U F (9), at ease. 

In the TFD formalism, all state vectors realized actually in the physical system have 
to be invariant under the tilde conjugation. Thus, not only the time-evolution but also 
all possible transitions of state vectors have to be invariant under the tilde conjugation. 

Requiring |a,7; 0) B given by Eq. (|36|) to be invariant under the tilde conjugation, we 
obtain 

|o;;0)b = |a, a;*;0) B 

= exp[aa f (0) + aa\6) - a*a(6) - a*a(6)} |O(0)) B . (39) 

From now on, we call this state the thermal coherent state [30], I3TJ [32]. Moreover, for 
simplicity, we assume a characterizing \a; 0) B to be always real. 

In the following paragraphs, we examine the physical meanings of the thermal vacua, 
|O(0)) B and |O(0)) F . 

First, we clarify the physical meanings of the thermal vacuum for the bosons |O(0)) B . 
We begin by considering the density operator defined on "H B , 

p B (0) = Tr^|O(0)) BB (O(0)|. (40) 

Then, we derive an explicit representation of p B (0) as follows: At first, from Eq. (|22|) and 
( 140|) . we obtain 

1 00 1 

PbW = — -2-E^B(6|5 n exp[(tanh0)a t a t ]|O,O) B 
cosh n\ 

x B (O,O|exp[(tanh0)aa](a t ) n |O) B . (41) 
Next, we apply the following relation to Eq. (I4ip : 

exp[-(tanh 0)a t a t ]5" exp[(tanh0)a t a t ] = [a + (tanh0)a t ] n . (42) 
Then using Eq. (lisp , we obtain 

1 oo 

= — r2^E tanh2n ^l^)BB(n| 
cosh n=o 

oo 

= (l-e-^^e-^I^BB^I. (43) 

n=0 

Looking at Eq. (1431) . we notice that p B (0) is an ensemble of quantum states {|n) B : 
n = 0,l,2,...}, into each of which n bosons (a-particles) are put. Moreover, the statistical 
probability of |n) B is given by (1 — e~^ e )e~ n/3t = Const, x e~^^ ne \ so that p B (0) represents 
a canonical ensemble of the Bose-Einstein distribution in thermal equilibrium. From 
these considerations, we understand that the a-particle represents the dynamical degree of 
freedom and the a-particle represents the thermal degree of freedom in |O(0)) B G H B £g>7i B . 
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Second, we clarify the physical meanings of the thermal vacuum for fermions \0(6)) F . 
We begin by considering the density operator defined on % F , 

p F (e)=Tr^\O(e)) FF (O(0)\. (44) 
From Eqs. f l30j) . fl34"j) and fj4~4"j) . we derive an explicit representation of p F (9) in the form, 
pp(6) = f(^|[cos^ + (sin^)c t £ t ]|0,0) FF (0,0|[cos^ + (sm6)cc]\n) F 

ne{0,l} 

= cos 2 0|O) FF (O| + sin 2 6>|1) FF {1| 

= r^^l°M0| + ^^|l>FF<l|. (45) 

Looking at Eq. (|45|) . we notice that p{6) F is an ensemble of quantum states {\n) F : 
n G {0, 1}}, into each of which n fermions (c-particles) are put. Moreover, the statistical 
probability of \n) F is given by (1 + e -^ e )- 1 e - n / 3e — Const, x e _/3 ( ne ), so that p{6) F repre- 
sents a canonical ensemble of the Fermi-Dirac distribution in thermal equilibrium. From 
these considerations, we understand that the c-particle represents the dynamical degree of 
freedom and the c-particle represents the thermal degree of freedom in |0(#))f G 'Hf£5>'Hf- 

For convenience of calculations that appear in the remains of this paper, using Eqs. ( TTBT) 
and ( 1T91) . we rewrite \a; 6) B given by Eq. ( )39|) as 

\a;9) B = exp[U B (6)(aa i + aa ] - aa - aa)Ul(6)]U B (6)\0,0} B 
= U B (8) exp[a(a^ — a)] exp[a(a* — a)] |0, 0)b 

= U B (9)\a) B \a) B , (46) 

where \a) B and |<5)b are the coherent states at zero temperature defined on T-L B and 7i B , 
respectively. 

Next, after the above preparations, we discuss how to construct the finite-temperature 
JCM according to the TFD formalism from the original JCM Hamiltonian. 

For example, we consider a system, which consists of bosons a and fermions c. We 
assume that the time-evolution of the system is reversible and it never causes dissipation. 
This implies that the Hamiltonian H, which is the Hermitian operator corresponding to 
the total energy of the system (the a-particles and the c-particles), is equivalent to the 
generator of the unitary operator for the time-evolution. 

Obeying the TFD formalism, we introduce the a-particle corresponding to the a- 
particle and the c-particle corresponding to the c-particle into the system. Because a 
and c are fictitious particles representing the thermal degree of freedom, the Schrodinger 
equation governing the time-evolution of the particles a and c never suffers form the 
thermal effects. This observation suggests that the Schrodinger equation for a and c 
never changes in spite of introducing a and c. From these considerations, we can conclude 
that the particles (a and c) are never coupled to the tilde particles (a and c) direct in the 
Hamiltonian. Hence, the total Hamiltonian based on the TFD formalism has to be a sum 
of the interaction terms of the particles (a and c) and the interaction terms of the tilde 
particles (a and c). 
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Here, we describe the total Hamiltonian for the TFD formalism as H . Then, letting 
the Schrodinger equation for the whole system consisting of the particles of a, a, c and c 
be invariant under the tilde conjugation, we can express H in the form, 

H — H — H, (47) 

where H = (H)~ [HI H21 H3J E]. If we give the Hamiltonian of the system H in the form 
of Eq. (14"T|) . the time-evolution of the whole system consisting of the particles and the 
tilde particles is reversible and it never causes dissipation. 

In Eq. ff4T|) . we emphasize the following: The Hamiltonian H is constructed from a 
and c. Then, because of H = (H)~ , the Hamiltonian H has to include both a and c. 
This implies that we have to introduce the temperature into both the a-particles and the 
c-particles. 

If the Hamiltonian of the total system is given by Eq. fT4Tj) . the thermal vacua |0(#))b 
given by Eq. (122]) and |0(6 i ))f given by Eq. fl34"|) are not dependent on time, so that they are 
stationary states. We can explain this fact as follows: Turning our eyes towards Eq. (1221) . 
we notice that the condensation of (aa)-pairs into |0(6 , )) B occurs. Thus, in the vacuum 
|0(#))b, the a-particle receives a phase factor exp[— i(H/h)t] and the a-particle receives 
a phase factor exp[i(H /%)t\. Then, these phase factors cancel out their effects with each 
other, and the (aa)-pair acts like a zero-energy boson. Hence, we understand that the 
thermal vacuum |0(6>))b is not dependent on time. We can apply a similar discussion to 
the thermal vacuum of the fermions, because we can observe the condensation of (cc)-pairs 
into |O(0)) F in Eq. (IMj). 

To obtain the JCM Hamiltonian H which takes the form given by Eq. (1471) . we rewrite 
the original Hamiltonian H of the JCM given by Eq. ([T|) as 

H = ^w (2c f c - 1) + hcoa^a + M(c ] a + ca ] ). (48) 

The Pauli matrices a z , a + and <r_ in the original Hamiltonian H given by Eq. ([1]) are 
replaced with the creation and annihilation operators of the fermions (2c^c — 1), c' and c 
in the rewritten Hamiltonian H given by Eq. (T481 . 

Extending the Hamiltonian H given by Eq. (|4"8"1) according to Eq. (|4T|) . we obtain 

H=(H + H 1 )-(H + H I ), (49) 

where 

H = —uo(2c^c — 1) + Tiuja)a } 

H\ = fiK(c^a + ca^), 

H = — w (2c^c — 1) + hcua^a, 

Hi = TiK{c ] a + ca ] ). (50) 
In this paper, we concentrate on examining the Hamiltonian H given by Eqs. (jUJ) and 

(EUD. 
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Here, we think about some Hamiltonians, which are proposed in the other works. 
Especially, we examine whether or not we can regard them as genuine JCM Hamiltonians 
based on the TFD formalism. 

Barnett, Knight and Azuma consider the following Hamiltonian in Refs. [To*! |2"2"]: 

h 

Hbka = —uq(J z + hoj(a'a — a)a) + hn,(a + a + cr_a'). (51) 

However, the Hamiltonian Hbka given by Eq. ( l5ip does not include the tilde operators 
corresponding to the atomic operators a z and a±, so that the states of the atom are 
always at zero temperature. Because the Hamiltonian Hbka given by Eq. (I5T]) does not 
introduce the temperature into the atom, we cannot regard it as a genuine Hamiltonian 
based on the TFD formalism. 

Fan and Lu propose the following Hamiltonian in Ref. 

ifpL = t^o^z + Tujj(cJa — a) a) 




/a — at r~ 77T f~, TT7 /at — a 
+nK[a + \ — r y a* a — a* a + y era — a*a\ (52) 



a — a 



The reason why Fan and Lu construct the Hamiltonian ifpL given by Eq. (|52|) is as follows: 
They find the commutation relations, 



[a) a — a^a, 





[a 1 a — a^a, \j — -} = ~\j — ~- (53) 

Thus, if we regard (a^a — a* a) as an extended number operator, we can think 
\J (at — a)/(a — dX) and \J{a — cfi) / (at — a) to be extended creation and annihilation op- 
erators, respectively. From these suggestions, Fan and Lu propose the Hamiltonian ifpL 
given in Eq. f )52|) . However, Fan and Lu's Hamiltonian ifpL does not include the tilde 
operators corresponding to the atomic operators a z and a±. Thus, it can not be regarded 
as a genuine Hamiltonian based on the TFD formalism. Moreover, in the Hamiltonian 
Hpi,, the operators a and are coupled direct to the tilde operators a and aL Thus, 
Fan and Lu's system suffers from dissipation during the time-evolution. To examine their 
system, we have to deal with non-equilibrium states. Because it is beyond the purpose of 
this paper, we do not involve ourselves in it. 

From now on, we examine the Hamiltonian H given by Eqs. (149p and ( )50l) . We can 
divide the Hamiltonian H into two parts as follows: 

H = h(C 1 + C 2 ), 

C\ = co[(c^c — c^c) + (a)a — a)a)\, 

C 2 = n^a + ccJ) - K^a + ca?) - Aw(c f c - c f c). (54) 

Then, we obtain a commutation relation [CijCy = 0. Moreover, we can diagonalize G\ 
at ease. Thus, we describe the time-evolution of the total system with the interaction 
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picture as follows: First, we write the state vector of the total system in the Schrddinger 
picture as |\l/g(i)). Second, assuming |^i(0)) = |\l/g(0)), we define the state vector of 
the total system in the interaction picture as = exp(iCit)\^s(t)). Hence, we can 

describe the time-evolution as \^i(t)) = U(t)\^i(0)), where U(t) = exp(—iC 2 t). 
The unitary operator for the time-evolution of is given by 

U(t) = exp[-iC 2 t] 
= U(t)®U(t) 

= exp [~ 2t ( ~tT Au/2 )] ® exp ^ ( ~tT Au/2 )]■ ^ 

In the right-hand side of Eq. (1531) . the first 2x2 matrix of the tensor product acts on Hp 
and the second 2x2 matrix of the tensor product acts on 1ri?. The 2x2 matrix U(t) 
appearing in Eq. (I55I) and the unitary operator for the time-evolution defined in Eq. 
are in the same form. Thus, the elements of the 2x2 matrices U(t) and U(t) appearing 
in Eq. (|55|) are given by Eq. (jSJ). 

4 Thermal effects of the period of the revival of the 
Rabi oscillations 

In this paper, putting the single cavity mode and the atom in the thermal coherent 
state | a; 9) B and the thermal vacuum |0(6 i ))f, respectively, at the time t = 0, we aim at 
examining the time-evolution of the JCM. We assume that the system consisting of the 
single cavity mode and the atom evolve in time without dissipation, and it maintains the 
constant temperature (3 all the time. As shown in Sec. [2, we can estimate that the period 
of the revival of the Rabi oscillations at zero temperature is around 27t|q!|/|k|. In this 
section, we discuss how the period changes at finite low temperature. We evaluate the 
thermal effects of the period in an intuitive manner. 

The parameter \a\, which characterizes the ordinary zero-temperature coherent state 
| a), is given by 

\a\ = ((alalia)) 1 / 2 . (56) 

Thus, we can guess that the parameter |a| varies with the thermal effects of the finite low 
temperature as 

\a\ {{a;9\a)a\a-6)) 1/2 . (57) 

On the other hand, using Eqs. (fT5l) . (|T6|) . (jTTj) . (|T9|) and (|46j) . and paying attention to 
a\a) = a\a) and a\ot) = a\a), we obtain 

(a; 9\a'a\a; 9) 
= (a\(a\U B (-9)<JaUl(-G)\a)\a) 
= (a|(a|a t (-^)a(-^)|a)|a) 

= a 2 e 29 + (l/4)(e 29 + e- 2e -2). (58) 
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From the above observations, we can expect the following: assuming a 2 3> 1 and 
9 <C 1, the thermal effects let the parameter characterizing the coherent state change 
under the low-temperature limit as 

\a\ -»■ ((a; 9\a^a\a; 9)) 1/2 ~ |a|e e . (59) 

Thus, we can expect that the period of the revival of the Rabi oscillations at finite low 
temperature varies as 

27r|or|/|K|-)>27r|a|e7l«|. (60) 

This phenomenon is confirmed by numerical calculations in Sec. [TUl 

The intuitive discussions given in this section is effective, when we can specify the 
whole system with the constant temperature (3. If the system is in a non-equilibrium 
state and the temperature (3 varies during its time-evolution, we cannot apply the above 
intuitive discussions to the system, so that Eq. fl60|) does not hold. 



5 The formulation of the perturbation theory 

In this section, we initially put the states of the atom and the cavity field in the thermal 
vacuum of the fermions \0(9)) F and the thermal coherent state \a; (3) B , respectively. Then, 
we formulate the time-evolution of the JCM based on the TFD discussed in Sec. [3] as the 
perturbation theory under the low-temperature limit. After formulating the perturbation 
theory here, we estimate the zero-th, first, second and third order corrections in Sees, d 
[HI [9] and Appendix |A3 To evaluate these correction terms, we make use of techniques for 
calculations developed in Ref. [22J. 

At first, we express the state of the system for t = in the form, 

|* I (O)) = |O(e)) P |a;0) B . (61) 

During the time-evolution of the system, we assume the atom and the cavity mode do 
not suffer dissipation and maintain the constant temperature /3. Thus, from Eqs. (fT8l) 
and (130]) . the parameters of the temperature for the fermionic atom Q(/3) and the bosonic 
cavity mode 6(f3) are given in the following forms, respectively: 

cos9(/3) = [1 + exp(-/3frw )r 1/2 , 

sin6(/3) = exp(-/3^ /2)[l + exp(-/3^ )]- 1/2 , 

cosh#(/3) = [1 - exp{-f3hu)]- 1/2 , 

sinhfl(/3) = [exp(phou) - 1]~ 1/2 , (62) 

where Uq represents the transition frequency of the two-level atom and u represents the 
frequency of the single cavity mode, as defined in Eq. ([1]). From Eqs. fH6|) and fl55|) . we 
obtain as 

\Mt)) = U(t)\%(0)) 

= [U(t)®U(t)}\0(Q)} F \a;9} B 

= [U(t)®U(t)][\0(e)) F U B (9)\a) B \a) B \. (63) 
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From Eq. f )34|) . we can rewrite |0(O))f in the form, 

|0(9)) F = cos6|0,0) F + sin6|l, 1) F . 



(64) 



Taking {\i,j)p : i,j G {1,0}} for the basis vectors of the four- dimensional Hilbert space 
%-p <8> "Hfj we can write down |0(6))p as a four-component vector, 



10(6))! 



( sin6 ^ 



V cos0 ) 



(65) 



where the components of the above vector are arranged in the order of |l,l)p, |1,0)f, 
|0,1)f and |0,6) F . 

Thus, writing |^i(t)) as the four- component vector, we obtain 



\Mt)) = [u{t)®u{t)\ 



( sinef> B W|a) B |a)B \ 



\^ cos6f/ B (6 | )|a)B|<5)B J 



(66) 



Moreover, expressing U(t) ® U{t) in the form of the 4x4 matrix, 

U(t) ® U(t) -- 



u n ll(t) u 10 U(t) 
u 01 U(t) u 00 U(t) 



U(t) 



(67) 



(6* 



where {u^ : i, j G {1,0}} and {uij : i, j G {1,0}} are given by Eq. ([5]), we can write down 
the four- component vector l^'i(t)) as the following explicit form: 



l*i(*)> 



^01 

V V'oo I 



(69) 



where 



^00 



(sin GmuMh + cos QuiqU W )U b 
(sin ©Mnttoi + cos QuioUoo)Ub 
(sin ©mqiMh + cos Ou 00 Ui )Ub 
(sin Qu 01 u 01 + cos Qu qUoo)Ub 



a) B \a) B , 
a) B |a) B , 
a) B |a) B , 
a) B \a) B . 



(70) 



Hence, the probability that we detect the ground state of the atom at zero temperature 
in the state of the total system |^i(t)) is given by 



P g (G,e;t) 



<(0,0|*i(0>|| 2 + ||f(0,1|* 1 (*))|| : 
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= nv>o6ir+iiv>oiir 

= B (a | B (at | U B (8) [ sin2 ©?4i M oiMoi «oi 

+ sin 6 cos e(i4iW 00 4i«oo + mS MoiMoo«oi) 
+ cos 2 Qul u OQ ul uo Q 
+ sin 2 Gmqi'Uoi^ii'Sii 

+ sin 6 cos e^MoottiiMio + mJ m im1 mh) 

+ cos 2 9Mj )0 MooSloMio]^B(6')|a)B|a)B- (71) 

Here, we pay attention to the following fact: Because U(t) given by Eq. fl68|) is a unitary 
matrix, we obtain 



UnUu + «01«01 


= 1, 


ttjoMio + MoO^OO 


= 1, 


tt^Uio + Moi^OO 


= o, 


«10«11 + Moo^Ol 


= 0. 


(I7T]) yields 





(72) 



P g (Q,9]t) = B (a\ B (a\U B (9)[sm 2 0mJiWoi + cos 2 0i4o«oo]%(#)l a )B|a)B 
= cos 2 0B(a|B(<S|^B(6')5 , i( a ^ a + c )^ r B(6 l ) |c^)b|ok)b 

+ sin 2 0B(a| B (a|^B(%2(aWc + i)U B (6)\a) B \a)B, (73) 



where 



gi{a)a + c) = Uq 



UoqUqo 



cos 2 



, r~ . . , sin 2 ( i/ata + c\n\t) 

(Vat a + c k t) + c ^-r L - L A 

v v 117 at a + c 



# 2 («^a + c+l) = i%«oi 

sin 2 (V at a + c\K,\t) ± 



a- 

a T a + c 



sin 2 (A/ata + c + ll/cli) r . + . 

= a f a + c+l -c. 74 

ata + c + 1 LV ' J v 7 

Using Eqs. (fTEj) and (IT7j) . we can rewrite P g (Q,9;t) given by Eqs. (1731 and fl?4|) as 

P 9 (9, 9; t) = cos 2 G(a| (a\ exp[6>(aa — a) aJ)\gi(w a + c) exp[— 9(aa — aJa/)}\a>}\a) 
+ sin 2 0(a|(a:| exp[6*(aa — a^at)]g 2 (a^a + c + 1) 

x exp[— 9(aa — aW)]|a)|<5:). (75) 
Thus, we obtain the perturbative expansion of P g (Q,9;t) in the small parameter 9{(i) as 

Pg (Q,9;t) = f: e -^P^(Q,9;t), (76) 

n=0 U - 
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where 

PW(9, 6- 1) = cos 2 0P#(t) + sin 2 0pW(t), (77) 



= 


{a 


/ ~ 1 / t i N 1 

{a\gi{a'a + c)\a)\a) } 


= 


(a 


(a\g2(a'a + c + 1) \a) a), 




(a 


(a [aa — (J a\ gi(a) a + c)] a) a), 




(a 


{& [aa — aW, ^(a a + c + 1)] a) <5), 


p${t) = 


(a 


(a [aa — aW, [aa — aW, ^i(a^a + c)]]|a)|a), 


pglw = 


(a 


(5 [aa — aW, [aa — a)a\ g 2 (a a + c + 1)]] |a) a), 



FjJ(t) = (a|(a| [aa — aW, [aa — a^, ai(a^a + c)]...] |a)|a), 



n-fold bracket 



PjjMt) = (a| (a | [aa — aW, [aa — aW, (^(a^a + c + 1)]...] |a) |a) 
, < 

n-fold bracket 

for n = 1,2,3,.... (79) 



Here, we pay attention to the following fact: The perturbative expansion given by 
Eqs. ( 176|) . ( 1771) . ( 1781) and ( 1791) is a power series in the small parameter 9 ({3). On the 
other hand, all the terms of the parameter 0(/3) included in the perturbative expansion, 
namely cos 2 G and sin 2 0, are expressed as explicit rigorous forms. Thus, we can strictly 
compute the functions of 0(/3) at ease in Eqs. (ITS]) and (!77j) . so that we do not need to 
worry about perturbative corrections of the parameter Q(j3). In this paper, we consider 
the power series in the small parameter 9(/3) to be the perturbative expansion under 
the low-temperature limit. In contrast, we do not regard 0(/3) as the parameter for the 
perturbation. 

Furthermore, the following trick lets actual computations of correction terms, that is 
to say, P^ n )(0, 9; t) for n = 1, 2, 3, be tractable. We can write down the functions gi(x) 
and g2(x) defined in Eq. ( 1741) as 

o , / — i i v sin 2 (v / xk|t) 
gAx) = cos 2 (v^kt) + c vv 1 1 \ 

X 

92(x) = W 1 1 \ x-c), (80) 

X 

so that we can rewrite each of them as the Taylor series at x = 0, 



g .(x) = 9i x ™ for ~oo<x< +00, (81) 

m=0 



where 

1 d m 

for i e {1,2}. (82) 



( m ) 

Qi — — 7~i — Qi[X 
yi m\dx m y ' 



x=0 
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Thus, we can rewrite Eq. f!79|) as 

CO 

Pgi{t) = X! 9i ( a l(^l — o^a\ [aa — aW, {a)a + c) m ]...] |a)|a), 

m=0 S ' 

n-fold bracket 

oo 

= X 92 — aW, [oo — aW, (a^a + c + l) m ]...] |ot)|a) 

— ' s / 

m=0 v 

n-fold bracket 

for n = 1,2,3,.... (83) 

In Sees. El [HI [9] and Appendix [A] using the perturbative expansion given by Eqs. (1761) . 
dZZD, (HE}, (Hi, flEQD, (EH), §2} and (JB3J, we compute P g (e,8;t). 

6 Comparison of the TFD formalism and the 
Liouville-von Neumann equation 

In the previous sections, we discuss a method for examining the time-evolution caused 
by the Hamiltonian H defined in Eqs. (|49[) and (15771) with the initial state |^(0)) = 
|O(0))f|«; 0)b given by Eq. (1511) according to the TFD formalism. This method is equiv- 
alent to solving the following Liouville-von Neumann equation: 

^p(t) = -^-[H,p(t)}, (84) 

where 



p(0) = Tr^[|*(0)><*(0)|] 

= Tr^[|0(6)) FF (0(e)|® \a;6) BB (a;6\], (85) 

and the Hamiltonian H appearing in Eq. (184"]) is given by Eqs. Jl]) and ( l4"8"|) . 

Both the Hamiltonian if based on the TFD formalism defined by Eqs. ( H9|) and (15 Op 
and the Liouville-von Neumann equation given by Eq. flSIJ) represent that the total system 
evolves in time with maintaining the constant temperature, so that it never suffers from 
dissipation and its time-evolution is reversible. Thus, we understand that the Hermitian 
operator corresponding to the energy of the total system is equivalent to the generator of 
the unitary operator for the time-evolution. 

Here, thinking about the Liouville-von Neumann equation given by Eq. (I84p . we divide 
H into the two parts C\ and C2 as shown in Eq. ([2J) and take the interaction picture. 
Assuming pi(0) = p(0), we introduce the density operator described in the interaction 
picture as 

Pl (t) = e lClt p(t)e- lClt . (86) 
Then, using the commutation relation [Ci,C 2 ] = 0, we obtain 

^(t) = -i[C 2 , Pl (t)}. (87) 
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From Eq. ( 187)) . we notice that we can rewrite p\(t) as 



Pl {t) = e-* C2 Vi(0y C24 



e- iC2t [|0(9)) FF (0(e)| ® |a; fl) BB (a; 0|]e iC2 ' 



Moreover, the probability that we detect the ground state of the atom at zero temperature 
is given by 

P g (Q,0;t) = F (0|Tr B [pi(t)]|0) F . (89) 

The physical meaning of Eqs. (188)) and (189)) is equivalent to the discussion developed from 
Eq. (JSU until Eq. flTHJ in Sec. 

Thus, comparing Eqs. (j75J and (189")) . we cannot find distinct differences between the 
TFD formalism and the Liouville-von Neumann equation. However, if we take the TFD 
formalism, we can express a physical quantity as a power series in 0(13) such as Eqs. (175)) . 
( 177)) . ( ITS)) and (179)) . Because of this advantage, the TFD formalism is superior than the 
Liouville-von Neumann equation for computing physical quantities actually. The reason 
why we take the TFD formalism in this paper for describing the JCM at finite temperature 
is the fact mentioned above. And this prescription is a new key point of this paper as 
compared with the other past works. 

In fact, if we rewrite Eqs. ( 188)) and ( 189)) as a low-temperature expansion without using 
the TFD formalism, we have to carry out the following calculations: 

Ti^[\a;0) BB (a;0\] 
= Tr n iU B (O)\a) B \a) BB (a\ B (a\U B (0)] 

= Tr^(exp[— (aa — a^)]\a} B \a} BB (a\ B (a\ exp[0(ad — a) a 1 )] 
= Tr^(|a)B|<3)BB(a;|B(<S| — 0[aa — era', \a) B \a) BB {a\ B {a\] 



2 

+ — [aa — a)a\ [aa — a)a\ |o;)b|5:)bb( c ' : |b(«|]] + •••)• (90) 

The above calculations are essentially equivalent to Eqs. ( ITS)) . ( 177j) . ( 1781) and ( 1791 . How- 
ever, the perturbation theory via the TFD formalism provides us a clearer insight and a 
more accurate understanding than the Liouville-von Neumann equation does. 

7 The zero-th order correction 

From Eqs. ( 174)) . (177)) and (1781) . we can write down the zero-th order correction as 

Pf)(6, 0- 1) = cos 2 QP${t) + sin 2 QP${t), (91) 

Pg?i(t) = (oi\gi(a^ a + c)\a) 

_ a 2 ~ a 2n r 2 , — — sin 2 (v / nTc|fi;|£) 1 

= e a > — cos 2 (Vn + c k \t) + c — — LJ -], 

^ n\ n + c 
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Pg$(t) = (a\g 2 (Ja + c + l)\a) 



n=0 Ul 



l2 ^^TM (n + 1)i (92) 



n=o nl n + c + 1 



Referring to Eq. (HDD, we note that P$(t) = P g (t) and ?f(0,0;t) = P g (t). 

For the convenience of calculations carried out in the remains of this paper, we define 
the following functions, each of which is represented as an infinite series: 



00 2n 

v 2 \ - a 



n=0 n! 



2 ~ « 2n r 9 . , -, ., sin 2 (A/n + c + Zkli). 

e"° V — cos 2 Vn + c + / k t) + c K — r ±A1 \ 

n\ 1 v 1 1 ' n + c + l 

.2 ~ a 2n 



n=o n - 



From the above definitions, we obtain the zero-th order correction terms as 

Pflit) = Q?(t), 

P${t) = Q?\t). (94) 

8 The first order correction 

From Eqs. (ITT)) and fl83|) . we can write down the first order correction as 

PW(6, 0; t) = cos 2 9PS(0 + sin 2 QP$(t), (95) 



00 

= E^ n) ( a l(«l[ aS - atst '( ata + c ) n ]l«)l«)' 

n=0 

00 

P i,2(*) = Y,92 n) (a\(a\[aa- a^a\ (a^a + c + l) n ]\a)\a) . (96) 



n=0 



From Eq. ( 1961) . we notice that we have to calculate the commutation relations, 

[ad- a)a\(a ( a + c) n ] for n = 0, 1, 2, .... (97) 
At first, we define the following three operators: 

A = a)($ — ad, B = a) a + c, C = a) a) + ad. (98) 
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[We pay attention to the fact that the operator C defined in Eq. (198]) is different from C%, 
C 2 , C\ and C 2 given by Eqs. © and (IMll.l Then, we obtain the commutation relations, 



= -C, 

[A, B 2 ] = A-2BC, 

[A,B 3 ] = -C + 3BA - 3B 2 C, 

[A,B 4 ] = A - ABC + 6B 2 A - 4B 3 C, 



(99) 



Next, we define the following two operators: 

ft = a* a', = ad, (100) 

and we obtain 

A = fi-i>, (7 = /* + />. (101) 
Using the operators fx and z>, we can rewrite Eq. (19"9~|) as the general form, 

[A, B n ] = (B - l) n jl - (B + l) n v - B n A for n = 1,2,3,.... (102) 

We can prove Eq. f |102[) with the mathematical induction as follows: First, we can confirm 
that Eq. ( 1102 j) holds for n = 1, at ease. Second, we assume Eq. ( 1102 j) holds for some 
unspecified number n(> 1). Third, we compute the commutation relation, 

[A,B n+1 ] = [A, B n ]B + B n [A, B] 

[{B - l) n fi -{B + \) n v - B n A)B + B n [A, B] 
= (B-l) n+1 fi-{B + l) n+1 i>-B n+1 A, (103) 

where we use fiB = [B — \)fi and vB = (B + 1)0. 
Thus, we obtain 

[ad — a^d\ (a)a + c) ra ] 
= -(a) a + c - l) n a)a ] + (a) a + c + l) n ad - (a) a + c) n (ad - a ] d)) 
= —a^d^(a^a + c) n + (a^a + c + l) n ad — (a^a + c) n ad 

+a ] d\a ] a + c+l) n . (104) 

Moreover, replacing c in Eq. f |104j) with (c + 1), we obtain the commutation relation, 

[ad — a^d\ (a) a + c + l) n ] 
= -a)d\a)a + c+ l) n + (a f a + c + 2) n aa 

-(a f a + c + l) n aa + a+a^a + c + 2) n . (105) 

Hence, using the relations a\a) = a\a) and d\a) = a\a), substitution of Eq. (j!04p into 
Eq. dnSD yields 

oo 

= E^! n Vl(a|[-a t a t (a t a + c) n + (a t a + c+l) n aa 

n=0 

22 



-(a)a + c) n aa + o)a\a)a + c + l) n ]|a)|a), 
(a|(a|[— (V($gi{a*a + c) + g\{wa + c + l)aa 
—gi(a)a + c)aa + cvcrgi(cra + c + l)]|a) |a), 
— 2a 2 (a;| (a| [(?i(a^a + c) — <7i(a'a + c + 1)] |a) |a) 



oo 2n 
„2 t — > tt 



= -2a 2 e- a2 J2-r[9i(n + c)- gi (n + c+l)] 

= -2a 2 [Qf\t) - it)}. (106) 
Similarly, substitution of Eq. fl 105[) into Eq. (!96|) yields 

P$(t) = -2a 2 e- a2 J2—[92(n + c+l)-g 2 (n + c + 2)} 

n=0 n - 

= -2a 2 [g( 0) (t)-Q( 1} (t)]. (107) 



9 The second order correction 

From Eq. ( 1831) . we can write down the second order terms as 

oo 

PgA^f) = g{ (at\(a\[aa — a)a\ \aa — a)a\ (a^a + c) n ]]\a)\a), 

n=0 

oo 

P i,2 } (*) = Y.9t\a\(a\[aa-a^\[aa-a ] h\(a ] a + c+ l) n ]]|a)|a). (108) 

n=0 

Looking at Eqs. (1981) and (11081) . we notice that we need to calculate the commutation 
relation [A, [A,B n ]]. 

From now on, referring to Eq. (11021) . we divide the commutation relation [A, [A,B n ]] 
into the following two parts and examine each of them: 

[A, [A, B n ]\ = R n + S n for n = 1,2,3,..., (109) 

where 

Rn = [A,(B-l) n ]fi-[A,(B + l) n ]v-[A,B n ]A, 

S n = (B-l) n [A,fi]-(B + l) n [A,u]. (110) 

According to Eq. (11091) , we divide one of the second order terms given by Eq. (I108P into 
two parts as 

oo oo 

P 9A(t) = ^g{ n \a\(a\R n \a)\a) + ^g{ n \a\(a\S n \a)\a). (Ill) 

n=0 n=0 

Here, we show P$(t) concrete example in Eq. (II lip . We understand obviously that 
we can compute PgJ(t) after the manner of Pg{(t). Thus, for simplicity, we concentrate 

(2) 

on evaluating Pg{(t) in the following paragraphs. 
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At first, we examine the part which includes {i? n } in Eq. f illip . From Eq. (11 101) . we 
obtain 



Ro 

Ri 
R2 
R 3 

R 6 



0. 

0, 

-2[A,B]C, 

-3[A, B 2 )C + 3[A, B)A, 

-4[A, B 3 ]C + 6[i, B 2 ]A - 4[A, B]C, 

-5[i, B 4 ]C + 10[i, B 3 ]A - 10[A, B 2 ]C + 5[A, B}A, 

-6[A, B 5 ]C + 15[i, 5 4 ]A - 20[i, 5 3 ]C7 + 15[i, 5 2 ]A - 6[i, B]C, 



(112) 



Thus, using Eq. fl 1 2 [) and the following formula: 



E 

n=0 



n + m + 1 \ (n+m+l) n+ i 



1 tf 



a; 



(m) 



, , m -</l(») -01 

ml ax m 
for m — 1,2, 3, ... 



(113) 



where gf 7 ^ is defined in Eqs. (IHTj) and (1821) . we can rewrite the part including {R n } in the 
second order term of Eq. (II lip as 



(a\(a\(-t( n V)9 ( r 2) [A,B^ 

n=0 V 1 / 



+ E 

n=0 
00 

-E 

n=0 



n t 3 U! n+3) [i,^]i 



n + 4 \ ( n+4 ) r 



_A^ + ijM-iF 3 C , + .--)l«>l«> 

00/ i \n 00 1 

J2 { -^F n fl-J2-F n u)\a)\a) 



n=0 



7V. 



n=0 



\(a\ (g x {B - l)/z - <&(B + l)z> - ffa (S)i)i|a)|a), 



where 



x=B-l 



x=B+l 



U / N 



A. 



(114) 



(115) 



In the derivation of Eq. (I114p . we use Eq. (11021) in an effective manner. The form of F n 
in Eq. fTTT5|) reflects Eq. ffT02|) . 

Using the operators e ±d ^ dx , we can rewrite the first term of the right-hand side of 
Eq. (JTUD as 

(a\(a\ {[e- d/dx gi (B - l)fi - e^^g^B + 1)0 - e~ d/dx gi{B)A]fi 
-[e d ' dx gi {B - l)jl - e d ' dx gi {B + l)u - e d l dx g x {B)A]P) \a)\a). (116) 
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Then, we apply the following technique to Eq. (j!16p : 

.=gi(x±i), (H7) 



n\ dx n 

where X is an arbitrary operator. Thus, we can rewrite Eq. (jllfip as 

(a\{&\ ([gi(B - 2)jj, - 9l (B)0 - 9l (B - 

-\gx(B)p. - gi{B + 2)0 - 9l {B + l)A]u) \a)\a). (118 



Next, we examine the part including {S n } in Eq. fillip . From Eqs. (HOOD and (110 II) . 
we obtain 

[i,A] = [i,*>] = LM = -A (H9) 

where 

D = a ] a + a ] a + 1. (120) 

In the remains of this section and Appendix |XJ we use Eq. fl 1 1 9 [) and the following 
commutation relations often without notice: 

[A, C] = -2D, [A, D] = -2(7. (121) 

Then, we can rewrite S n given by Eq. ( 11101) as 

S n = -((B-l) n -(B + l) n )b. (122) 



[In Eq. (57) of Ref. |22j, a calculation concerning S n is wrong.] Thus, we can write down 
the part including {S n } in the second order term given by Eq. (Hill) as 

- (a\(a\U(B - 1) - 9l (B + l))D\a)\a). (123) 



Putting together Eqs. fillip . (j!14j) . (11181) and (11231) . we can write down the whole of 
the second order term as 

= (a\(a\[ gi (B - 2)fi 2 + 9l (B - 1)(-D -jiA- Afi) 
+ gi (B)(A 2 -vpL- piv) + g x {B + 1)(D + vA + Av) 
+ 9l (B + 2)v 2 ]\a)\&). (124) 

Here, to compute Pg 2 } (t) given by Eq. (I124p . we arrange jl in the left side of the product 
of operators and v in the right side of the product of operators. For the arrangement of 
operators, we carry out the calculations, 

{B-2) n fi 2 = jl 2 B n , 
(B-l) n (-D-fiA-Afi) = -2fi 2 {B + l) n + 2fiB n i), 

B n (A 2 -vpL- jlv) = fi 2 {B + 2) n - 4£i(B + \) n v - 2B n D + B n 2 , 
(B + l) n (D + OA + AO) = 2p,{B + 2) n - 2(B + l) n 2 + 2{B + l) n D 

for n = 1,2,3,.... (125) 
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Substitution of the above relations into Eq. (11241) yields 

P$(t) = (a\(&\[^g 1 (B)-2fi 2 g l (B + l)+2jlg 1 (B)v + fi 2 g 1 (B + 2) 
-4m(B + l)v - 2 9l (B)D + 9l (B)u 2 + 2ft 9l (B + 2)z> 
-2^(5 + l)i> 2 + 2^(5 + 1)1) + 9l (B + 2)i> 2 ]|a)|a> 
= Aa 4 (a\(a\[ gi (B) - 2 9l (B + 1) + ^(B + 2)}\a)\a) 

-2(a|(a|</i(S)D|a)|a) + 2(a|(% 1 ( J B + l)D|a)|a). (126) 

Moreover, preparing the following formula: 

(a\(a\ gi (B)D\a)\a) 

oo oo Q,n+m 

= (a|(a|e~ Q2 ^ ^ -==(n + m + l)gi(n + c)|n)|m> 
n=om=o vn!m! 

oo (^2n oo Q,2n 

= e"" 2 ^ — ^i(n + c) + e~° 2 ^ — n#i(n + c) 

oo oo 2(n+m) 

oo ry 2n oo -,271 

= (l + a 2 )e - Q2 E^i(^ + c ) + « 2e " a2 E^i(^ + c + 1 )' (127) 
we arrive at the final representation of Pg{(t) as 

oo ^ ,2n 
,(2)/.s ■ * 2 ^ OL 



PiJ(t) = 4a*e- a J2 —Mn + c) - 2 9l (n + c + I) + 9l (n + c + 2)} 

n=0 n - 

o 00 a 2n o 00 r> 2n 

-2(1 + a 2 )e- a2 E ^r<?i(™ + c) - 2« 2 e" a2 £ ^r<?i(n + c + 1) 



n=0 H[ nT0 n! 

00 a 2n , , , „2 ^ a 2 " 



+2(1 + a 2 )e~ a V -r^i(n + c + 1) + 2a 2 e - Q V ^i(n + c + 2) 

= 2(2a 2 + l)(a + l)(a-l)Q { °\t) 
-2(2a 2 + l)(2a 2 -l)QS 1) W 

+2a 2 (2a 2 + l)Qf ) (t). (128) 

Similarly, we obtain P^it) as 

P$(t) = 2(2a 2 + l)(a + l)(a-l)Qf\t) 
-2(2a 2 + l)(2a 2 -l)Q^\t) 

+2a 2 (2a 2 + l)Q ( 2 } (t). (129) 

10 The numerical calculations 

In this section, we show numerical results for the atomic population inversion obtained 
with the third order perturbation theory under the low-temperature limit. In Sees. [7J |BJ 
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Eland Appendix EJ we obtain {Pj${t), P$(t) : n G {0, 1, 2, 3}} in the form of Eqs. AMI). 
( 1T06|) . ( llOTj) . f lT28|) . ( 11291) . f Tl69j) and ([110). Thus, from Eqs. ©, (EBJ and (EZD, we can 
calculate (cr z (t)) as the third order perturbation theory, 

<*,(*)) = 1 - 2[cos 2 009) £ ^f-P^(t) + sin 2 9(/3) £ (130) 

Figured] shows the atomic population inversion (c z (t)) given by Eq. ( 1130 j) as a function 
of the time t with a — 4, c — 1, k — 1 and 0(/3) = #(/3) = 0. Figure [3] shows the atomic 
population inversion (c z (t)) given by Eq. ( I130p as a function of the time t with a = 8, 
c = 1, k — 1 and Q(/3) = 9(0) = 0. Carrying out numerical calculations for Figs. [T] 
and El we replace the summation J2%Lo i n Qi(t) with 5^= , so that we compute the 
sum of first one hundred and one terms in the series. (In this section, whenever we carry 
out numerical calculations of Q%(t) and Q% (t), we replace their summation J2^Lo with 
J2n=o-) I n Fi§ s - E an d El we assume the system to be at zero temperature. Thus, the 
graphs in Figs. [T] and [3] do not suffer from thermal effects. We can observe the collapse 
and the revival of the Rabi oscillations obviously in these graphs. 

Figure HI shows the atomic population inversion (a z (t)) given by Eq. ( I13(jp as a function 
of the time t with a = 4, c = 1, k = 1, 9(f3) = vr/32, u = 2 and u = 4. From Eq. (152]) . 
we obtain 

= arctanh(e- /3?iw/2 ), 
0(/3) = arctan(e- /3 ^ o/2 ), (131) 

so that the relation exp(— 2/3h) = tanh[6 l (/3)] = tanh(7r/32) holds. Thus, we can derive 
the following relation: 

Q((3) = arctan(e-^) = arctan[tanh 1/2 (^(/3))] = arctan[tanh 1/2 (vr/32)]. (132) 

Because the system of Fig. [2] evolves in time with maintaining constant low temperature, 
its time-evolution is under the thermal effects. Comparing the graphs shown in Figs. [H 
and [21 we notice that the period of Fig. [2] is longer than the period of Fig. [Q Thus, we 
can suppose that the thermal effects let the period of the revival of the Rabi oscillations 
become longer. 

Figured] shows the atomic population inversion (a z (t)) given by Eq. f ll30p as a function 
of the time t with a = 8, c = 1, k = 1, 9((3) = 7r/60, uq = 2 and u = 4. Then, in a similar 
manner for obtaining Eqs. flTBTD and fTT32l . we achieve 6(/3) = arctan[tanh 1/2 (#(/3))] = 
arctan[tanh 1//2 (7r/60)]. Comparing the graphs shown in Figs. [3] and HI we notice that the 
period of Fig. @] is longer than the period of Fig. [31 so that we can suppose that the 
thermal effects let the period of the revival of the Rabi oscillations become longer. 

When we take a = 4, c = 1, K — 1, < t < 20tt and 9(0) = tt/32, a numerical value 
of each order perturbation correction varies as shown in Table [H On the other hand, 
when we take a = 8, c = 1, k — 1, < t < 40ir and 9(0) = 7r/60, a numerical value of 
each order perturbation correction varies as shown in Table [2] Turning our eyes towards 
Table [H we observe that the contribution of the third order correction is nearly equal to 
a half of the contribution of the second order correction in the perturbative expansion. 
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Figure 1: The atomic population inversion (cr z (t)) as a function of the time t obtained 
from numerical calculations of Eq. ( J13Uj) with a = 4, c = 1, n = 1 and ®(/3) = 9(0) = 0. 
Looking at the graph, we estimate the time scale of the initial collapse and the period of 
the revival of the Rabi oscillations at unity and 8n around, respectively. 



1 r 




-1 1 1 1 1 ' t 

5tt IOtt 15 77 20tt 



Figure 2: The atomic population inversion (cr z (t)) as a function of the time t obtained 
from numerical calculations of Eq. ( I130P with a — 4, c = 1, k — 1, Uq — 2, uj — 4, 
0(/3) = 7r/32 and Q(/3) = arctan[tanh 1//2 (7r/32)]. Looking at the graph, we estimate the 
time scale of the initial collapse and the period of the revival of the Rabi oscillations at 
unity and 87re n ^ 32 ~ (8.83)7r around, respectively. 
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Figure 3: The atomic population inversion (cr z (t)) as a function of the time t obtained 
from numerical calculations of Eq. ( J13Uj) with a = 8, c = 1, k = 1 and 0(/3) = = 0. 
Looking at the graph, we estimate the time scale of the initial collapse and the period of 
the revival of the Rabi oscillations at unity and 16n around, respectively. 
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Figure 4: The atomic population inversion (cr z (t)) as a function of the time t obtained 
from numerical calculations of Eq. (I130p with a — 8, c = 1, k — 1, Uq — 2, uj — 4, 
0(/3) = 7r/60 and Q(/3) = arctan[tanh 1//2 (7r/60)]. Looking at the graph, we estimate the 
time scale of the initial collapse and the period of the revival of the Rabi oscillations at 
unity and Wne^^ 60 ~ (16.9)7T around, respectively. 
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Table 1: The ranges of numerical values of the perturbation corrections with a — 4, c = 1, 
k — 1, < t < 20-7T and 0(j3) = 7r/32. The estimations of the minimum and the maximum 
in every row of the table are based on values of each correction term, which we obtain 
numerically at equally spaced intervals At = 20n x 10~ 4 during < t < 20n. 



correction term 


min 


max 


(l/2)0(/3) 2 P$(O 
(1/2)6103)^2(0 
(l/6)6(f})*P$(t) 

{i/sw?p?kt) 


-0.227 

-0.208 

-0.108 

-0.0996 

-0.0441 

-0.0478 


0.199 

0.225 

0.102 

0.109 

0.0483 

0.0467 



Table 2: The ranges of numerical values of the perturbation corrections with a = 8, c = 1, 
k — 1, < t < 407T and 9((3) = n/60. The estimations of the minimum and the maximum 
in every row of the table are based on values of each correction term, which we obtain 
numerically at equally spaced intervals At = 407r x 10~ 4 during < t < 40-7T. 



correction term 


min 


max 


0W)P^(t) 

mp$(t) 

(l/2)Wp£?(t) 
(l/2)<?(/3) 2 P$(f) 
(1/6)^/3)^(0 


-0.246 

-0.248 

-0.125 

-0.127 

-0.0566 

-0.0565 


0.247 

0.245 

0.128 

0.126 

0.0570 

0.0567 
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ln(T) 
3.36 



3.32 



3.28 



3.24 

tt/128 n/64 3^/128 tt/32 

Figure 5: The period of the revival of the Rabi oscillations T(9) plotted as a function of the 
parameter of the temperature 9((3). The points are obtained from numerical calculations 
of the third order perturbation theory with taking a = 4, c = 1, k = 1, cuq = 2, u = 4 and 
< 9(0) < vr/32. In the graph, the vertical axis is scaled logarithmically as ln[T(6 1 )] and 
the horizontal axis is scaled linearly as 9(0). Fitting the points with the linear function 
according to the least-squares method, we obtain ln[T(#)] = 3.25 + (0. 988)0, which is 
drawn in the graph. 



From Table [T], we consider the perturbative expansion to be reliable for 9(0) = 7r/32. 
Thus, taking a = 4, c = 1 and k = 1, we can conclude that the third order perturbation 
theory is effective for the parameter < 9(0) < n/32. We notice that a similar thing 
happens in Table [2j as well. Thus, taking a = 8, c = 1 and k = 1, we can conclude that 
the third order perturbation theory is effective for the parameter < 9(0) < n/60. 

In Figs. El and El we plot the period of the revival of the Rabi oscillations T(9) as 
a function of the parameter of the temperature 9(0). The points in Fig. El are obtained 
from numerical calculations of the third order perturbation theory with taking a = 4, 
c = 1, k — 1, ojq — 2, lu — 4 and < 9(0) < n/32. The points in Fig. El are obtained 
similarly with taking a — 8, c = 1, k — 1, ojq — 2, oj — 4 and < 9(/3) < 7r/60. In 
Fig. EJ we compute the period T(9) numerically as follows: First, we calculate (cr z (t)) 
given by Eq. (11301) for a certain 9(j3). For every point of Fig. El taking the interval of 
the time At = 5tt x 10~ 4 , we obtain (u z (t)) at each time step during 157r/2 < t < IOtc. 
We write the time at which (a z (t)) takes the maximum value as t max and write the time 
at which (cr z (t)) takes the minimum value as t m i n . Second, we obtain the period T(9) 
with taking T = (t max + t min )/2. [For example, taking 9 = n/32, we obtain £ max ss 28.52 
for (cr 2 (t max )) 0.3923 and t min ks 28.86 for (er 2 (t min )) ps —0.4763. Thus, we obtain 
T(7r/32) ~ 28.69.] The points of T(9) in Fig. El are obtained in a similar manner with 
taking the interval of the time At = IOti x 10~ 4 and carrying out calculations of (a z (t)) 
at each time step during 157r < t < 207T. 

In the graphs of Figs. El and El the vertical axes are scaled logarithmically as ln[T(0)] 
and the horizontal axes are scaled linearly as 9(0). In the graph of Fig. El the points 
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Figure 6: The period of the revival of the Rabi oscillations T(6) plotted as a function of the 
parameter of the temperature 0({3). The points are obtained from numerical calculations 
of the third order perturbation theory with taking a = 8, c = 1, k = 1, co = 2, cj = 4 and 
< 0(0) < 7r/60. In the graph, the vertical axis is scaled logarithmically as ln[T(6 1 )] and 
the horizontal axis is scaled linearly as 0(0). Fitting the points with the linear function 
according to the least-squares method, we obtain ln[T(6 1 )] = 3.92 + (1.07)0, which is drawn 
in the graph. 



form groups consisting of twos, threes and fours, so that they appear in the shape of the 
stairs as ln[T(0)] increases gradually. The reason why the points appear in the shape of 
the stairs is as follows: The atomic population inversion (a z (t)) is a bunch of the Rabi 
oscillations whose period is 7r/(|a||/c|) — vr/4 around. (We obtain this approximation in 
Sec. (21) At the same time, it shows the revival of the amplitude envelope with the period 
T(0) ~ 27r | ct | e 6 */ 1 | = 8ne 8 around. Thus, calculating t max and t m i n numerically, the rapid 
Rabi oscillations give us the smallest interval measurable as about tt/8, which is the half 
of the period of the Rabi oscillations. This resolution of the time lets the points in Fig. [5] 
form the shape of the stairs. 

Contrastingly, the points in Fig. [6] do not appear in the distinct shape of the stairs. 
This is because the resolution of Fig. [HI is finer than that of Fig. [5J Indeed, in Fig. [HI the 
period of the Rabi oscillations is given by tt/8 around, so that the resolution of T(0) is 
nearly equal to tt/16. 

Fitting the points in Fig. [5] with the linear function according to the least-squares 
method, we obtain 

In [T(0)] = 3.25 + (0.988)0. (133) 

We can interpret the above result as 

T(0) ~ e 3 - 25 x e^ e ~ 8ne 6 , (134) 

which reminds us of Eq. (!60l) . On the other hand, fitting the points in Fig. [6] with the 
linear function according to the least-squares method, we obtain 

In [T(0)\ = 3.92 + (1.07)0. (135) 
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We can interpret the above result as 

T{6) ~ e 3 - 92 x e( L07 ) e ~ 16vre e , (136) 
which also reminds us of Eq. (!60|) . 

11 Thermal effects of the counter-rotating terms 

In this section, we address thermal effects of the counter-rotating terms. Because this 
topic is difficult and includes subtle problems, we treat it with an intuitive manner. 

First of all, we have to go back to the derivation of the JCM. At the beginning, we 
consider the Hamiltonian for a magnetic dipole in a magnetic field, and we obtain 

h 

H = —u)q(J z + hcoa^a + hn(a^ + o~ + )(a + a'). (137) 

Assuming near resonance u ~ ujq, the interaction terms a + a and cr_at are practically 
independent of the time t, while the terms <r_a and <r+at vary rapidly at frequencies 
±(a>o + w)- Then, applying the rotating wave approximation to Eq. (1137P and removing 
the term hn(a + a^ + cr_a), we obtain the Hamiltonian of the JCM written down as Eq. (P). 

As mentioned above, the rotating wave approximation is used often in the field of 
the quantum optics. However, it is shown that the rotating wave approximation cannot 
always be a good treatment, and sometimes it causes serious defects. Ford et al. examine 
the Hamiltonian for an oscillator of the frequency ujq interacting with a reservoir and its 
rotating wave approximation [3"3"| 151] . The Hamiltonian of the original model is given by 

H = huoata + £ + ( a + flt ) £ x A b j + & ])> ( 138 ) 

3 3 

where [a, at] = 1, [bj, 6]] = 1 Vj, and its rotating wave approximation is given by 

-f^RWA = huoa^a + ^2 hujjbjbj + \j(aJbj + ao]). (139) 
3 j 

Then, the Hamiltonian -Hrwa defined in Eq. (I139P causes the following problem: The 
expectation value (the energy) of -Hrwa has no lower bound, so that we cannot specify the 
ground state. Thus, we have to think the system described with -£Trwa to be unphysical. 

As explained above, the rotating wave approximation sometimes manifests anomalous 
aspects. Someone might complain that the rotating wave approximation brings us the 
JCM that is an exactly soluble quantum mechanical model for arbitrary Aw and k. How- 
ever, the JCM also has a defect, which we cannot neglect. Here, we think around the 
eigenvalues and the eigenvectors of the Hamiltonian of the JCM given by Eq. (JT]). They 
are written down as follows [3]: 

E n> i = h[u(n+ ^) + A n ], 

En t2 = h[u{n + -) - A„] for n = 0, 1,2, 

£o,o = -~w , (140) 



33 




-28 



Figure 7: E 0)0 and {E n ^ '■ n — 0, 1, 2, 24} as functions of k, where < k < 10. Plotting 
them as graphs, we assume ujq = u = 1 and % — 1. Because never can be the 

ground-state energy, we do not plot them in this figure. Looking these graphs, we notice 
the following facts. When k = 0, the ground-state energy is equal to E^o- Ori the other 
hand, when n = 10, the ground-state energy is given by £'24,2- In fact, these graphs show 
that the ground-state energy changes from E 0t0 to E n>2 for n ^> 1 gradually as k becomes 
larger. 



\tp(n, 1)) = cos0 n |n + 1) P |#) A + sm0 n |n) P |e) A , 

\<f(n, 2)) = — sin^Jri + 1)p|(7)a + cos6' n |n)p|e)A for n = 0, 1, 2, 

|o,o) = |o>p|0>a, (141; 



where 



Au = lo — wo, and 



A n = ^)2 + ^(n + l), (142) 

tan^= "^+ r -. (143) 
(Aw/2) + A n 

Looking at Eq. fll40p . we notice that the ground state changes from |0, 0) to \ip(n, 2)) 
for n ^> 1 gradually as \k\ becomes larger. To confirm it numerically, we plot _Eo,o and 
{E n ^ '■ n = 0, 1, 2, 24} as functions of k in Fig. [71 At the same time, an excitation 
energy, which is required to promote the JCM system from the ground state to the first 
excited state, becomes smaller rapidly as \k\ — > 00. To confirm it numerically, we plot the 
excitation energy as a function of k in Figs. [S] and O 

From the analyses performed in Figs. [7J [S] and O we can conclude as follows: If we 
take a large value of the ground state of the JCM contains many photons. Then, the 
excitation energy takes a small value. These properties of the ground state of the JCM 
relate to the uncertainty principle AiVA0 > (1/2). Because AE decreases exponentially 
as \k\ becomes larger as shown in Figs. El and El the system of the JCM is able to jump 
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Figure 8: An excitation energy AE, which is required to promote the JCM system from 
the ground state to the first excited state, as an function of k, where < k < 10. Looking 
at this graph, we notice that AE oscillates and its amplitude decreases rapidly. When 
AE = 0, the JCM system has the degenerate ground states. 




Figure 9: An excitation energy AE, which is required to promote the JCM system from 
the ground state to the first excited state, as an function of k, where < k < 10. The 
graph uses the logarithmic scale on the vertical axis and the linear scale on the horizontal 
axis. In Fig. [HJ we show that the system has the degenerate ground states at certain 
values of the parameter k. Because In AE — > — oo as AE — > 0, we cannot plot small AE, 
which is nearly equal to zero, in the graph. Looking at this graph, we notice that the 
amplitude of AE, which oscillates in the parameter k, decreases exponentially. 
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from the ground state to excited states at ease for \k\ ^> 1. Thus, the fluctuation of 
the number of photons AN becomes very larger. Hence, according to the uncertainty 
principle ANAcp > (1/2), the system of the JCM around the ground state acquires very 
small fluctuation of the phase of each photon, so that A0 — > 0. 

However, the ground state that contains a large number of photons with small fluctua- 
tion of the phase seems not to be practical. We may realize the ground state for AiV 3> 1 
and A(p ~ in the laboratory by using a two-level atom in the cavity field, which is 
induced by a very strong laser beam. 

From the viewpoint explained above, we cannot regard the JCM derived with the 
rotating wave approximation as a proper model in the field of the quantum optics. Hence, 
the JCM is valid and has physical meanings if and only if a near resonance u ~ uq is 
assumed and \k\ is small enough. 

To overcome the defects of the rotating wave approximation in the JCM, some re- 
searchers try to extend and generalize the JCM. Ng et al. investigate the two-photon 
JCM and the intensity-dependent JCM with the counter-rotating terms [3~5| 136]. In these 
models, the nonlinearity of the interaction between the two-level atom and the cavity field 
is emphasized. 

In general, it is very difficult and complicated to evaluate the contributions of the 
counter-rotating terms in the JCM. Feranchuk et al. study the Schrodinger equation, 
whose Hamiltonian is given by Eq. ( I137p . numerically |37j . 

Especially, Phoenix presents several perturbative approaches to investigate this prob- 
lem. Here, we review one of his perturbation methods, which is called short time expansion 
of the inversion. First, we begin with the Heisenberg picture of a z , 

a z {t) = eMj-Ht)a z exp(-~#t), (144) 

where the Hamiltonian is given by Eq. ( I137p . Moreover, we assume u = Uq, so that we 
consider the optical resonance. 

Second, we expand Eq. ( I144p in a power series in t and neglect third-order terms [58] . 
So that, we obtain 

n 2n n 

= l-2(/tt) 2 (a + a t ) 2 + C(t 3 ). (145) 

Third, we assume the initial state as |o;)p|e)A, where a = y/Ee 1 ^, and substitute it into 
Eq. (I145p . Finally, we obtain 

(cr z {t)) a = A(e\p(a\a z (t)\a) P \e) A 

= 1 -2(/tt) 2 (4ncos 2 0+ 1) + C(t 3 ). (146) 

Because Eq. ( 1146P is valid for < t <C 1, we can expect it to describe the initial collapse 
of the Rabi oscillations. 

We note that Eq. (I146p depends on the phase 0. This characteristic can always be 
found in any perturbative expansion of (<J z (t)) a . (This fact is indicated by Phoenix first.) 
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In Sec. HI we give the intuitive discussions about the thermal effects of the JCM, and 
we obtain Eq. fl59l) under the low-temperature limit. Hence, we can add the thermal 
effects to Eq. jUB} as 

M*))a = 1 - 2(Kt) 2 (4ne 20 cos 2 + 1) + 0(t 3 ), (147) 

where 9 is given by Eq. f JT8|) . To examine whether or not Eq. (j!47p holds remains to be 
solved in the future. 

In this section, we argue only the short time expansion of (a z (t)) a with the counter- 
rotating terms. We point out that to examine long time behaviour of the JCM with 
counter- rotating terms is very difficult even if we use perturbative techniques. 

12 Discussion 

Turning our eyes towards the graphs shown in Figs. and El we observe that the period 
of the revival of the Rabi oscillations becomes longer as the temperature rises. This 
phenomenon is predicted form an intuitive discussion in Sec. HJ In Sec. [TUJ we confirm 
this phenomenon (or this expectation) with numerical calculations based on the third 
order low-temperature expansion. 

Why is the low-temperature expansion in 6((3) given by Sec. effective for a pertur- 
bation theory? The reason why is as follows: The thermal coherent state is defined in 
Eqs. (139]) and (146]) . This definition is suitable for the low-temperature expansion because 
of the Baker-Hausdorff theorem [3]. 

A The third order correction 

In this section, we give details of calculations of PgJ{t) and PgJ(t) defined in Eq. (183]) . 
From Eq. (153"]) . we can write down the third order correction terms as 

oo 
n=0 

[ad — a* a', (a/ a + c) n ]]]|a)|<5), 

oo 

Pgiif) = ^2 fi , 2™' ) ( a l( ( 5|[ a <2 — a^a\ [ad — a) a), 

n=0 

[aa-a + a t ,(a t a + c+ l) n ]]]\a)\&). (148) 

As mentioned in Sees. [7J[H]and|Hl w e can obtain pfty(t) after the manner of P$(t). Thus, 

from now on, we compute P$(t). For the convenience of calculations carried out in the 

remains of this section, according to Eq. f)109p . we divide Pg 3 i(t) into the following two 
parts: 

oo 

p ®(t) = -T,9{ n \a\(a\[A[A[AB n ]]}\a}\a) 

n=0 

oo oo 

= -E^ n ^«K 5 l[A^]|«}|«}-E^ n) («K 5 l[A^]i«}|a). (149) 

n=0 n=0 



37 



First, we examine the part which includes {[A, Rn)} in Eq. (I149p . From Eqs. (IllOp 
and (I112p . after slightly tough calculations, we obtain 

[A,R ] = 0, 
[A, Ri] = 0, 

[A,R 2 ] = -2(R 1 + S 1 )C-2[A,B][A,C], 
[A,R 3 ] = -3( J R 2 + 5 2 )C'-3[i,5 2 ][i,C , ] + 3(i?i + 5 1 )i, 
[A,R 4 ] = -A(R 3 + S 3 )C-A[A,B 3 }[A,C} + 6(R 2 + S 2 )A-A(R 1 + S 1 )C 
-A[A,B][A,C], 

[A,R 5 ] = -5(R i + S 4 )C-5[A,B 4 ][A,C] + 10(R 3 + S 3 )A-10(R 2 + S 2 )C 

-10[i, J B 2 ][i,C'] + 5( J R 1 + 5 1 )i, 
[A, Re] = -6(R 5 + S 5 )C-6[A,B 5 ][A,C} + 15(R 4 + S 4 )A-20(R 3 + S 3 )C 

-20[i, B 3 ][A, 6} + 15(4 + S 2 )A - 6(4 + St)C - 6[A, B] [A, C], 

(150) 

Thus, we can rewrite the part including {[A, R n ]} in Eq. (I149p as 

oo 00 / i 9 \ 

J2g{ n) [AR n ] = -E )9 { " +2) (Rn + i + S n+1 )C 



n=0 n=0 



-El n ^ 2 )g[ n+2) [A,B^}[AC] 



n=0 

oo 



+ E( n t 3 )g{ n+3 \Rn + i + s n+1 )A 



n=0 

oo 



-E( n 3 4 )^ +4) (^ + i + ^ + i)^ 

-£( n tM^ n+4) [i,5" +1 ][i,(7] 



n=0 

oo 



+ E( n t 5 U! n+5) (4 + i + 4 + i)i 



n=0 

oo 



E( "t 6 K +6) (4 +1 + 4 +1 )c 



n=0 

oo 



-E] n 5 6 )r(i + i + U[Aq 

-.".T° (151) 

In the following paragraphs, we examine the terms in the right-hand side of Eq. (11511) . 
one by one. 

Here, we think about the first term of the right-hand side of Eq. (I15ip . At first, 
referring to Eq. (I112p . we calculate the part including {R n+ i} in the first term of the 
right-hand side of Eq. (11511) as 

n=0 V / 
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-E 

n=0 

oo 

+ E 



x ' n + 2 
1 



n=0 

oo 



-E 

n=0 
oo 

+ E 



n=0 



n + 3 
2 

n + 4 
3 

72 + 5 

4 



(' n + 4 ^f +4 )[i, J B"+ 1 ]i 

^ + 5 \ g (n+5)^,B n+1 ]C 

( n + 6 )g[ n+6) [A,B^)A 



We prepare the following formula: 

n + m + 1 



E 

n=0 \ 
1 (f" 



r?) 



U + m + 2 \ ( n+rn+ 2) n+l 

I I yi x 



^(x) — (m + l)^ m+1 ^ for m — 1, 2, 3, .. 



ml dx m ^ 

From Eq. (11021) and the above formula, we can rewrite Eq. (11521) as 



n=0 



. P---j-sfi{x) 
x=b-i ax 

, 1 d 2 



x=b+i ax 



x=B-l 



3! dx 3 ' 

+... 

(e" d/ ^( 
+ (- e '/«Vi(a 
+ (-e- d ^g[(x) 



x=B-l 



x=B-l 



2! dx 2 

3! dx 3 
1 rf 4 
ttdx 1 



A)C 

x=B ' 



x =b+i 2! dx 2 
1 d 3 , 



x=B 



A) A 



x =b+i 3! dx 3 



g'A< 



x=B 



A)C 



x=B+l 
—d/dXgt 

+ g' 1 (B-l)+e d > d *g[(x 



4! dx 4 ' 



x=B 



A) A 



_^ (jB _l)_ e -^ (x) +g[(B))jl 2 
x=B-l x=B J 



+ (e d/dx g[(x 



x=B-l x=B 

+ g [(B + l) + e- d ^g[(x 

x=B+l 

-g'^B + V-e^g'^x 



-g'x{B))tiv 



x=B 



g[(B)m 



x=B+l 



x=B 



+ g'i(B))u 



g , 1 (B-2)-2g[(B-l)+g' 1 (B))fi 2 
+ (g[{B - 1) - 2g[(B) + g[(B + l))fw 
+ (g[(B - 1) - 2g[(B) + g[(B + l))vft 
+ (g[(B)-2g' 1 (B + l)+g[(B + 2 



v 2 . 
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Next, from Eqs. (1113)) and (j!22p . we can compute the part including {S n +i} in the 
first term of the right-hand side of Eq. (j!5ip as 



n=0 

OO 



-E( n ! 2 ui n+2) ((s-ir +i -(s+ir +i )i) 

n=0 V / 

-fe(B-l) -^(fl + !))£>. (155) 



From Eqs. ( HU2p and (11 13)) . we calculate the second term in the right-hand-side of 
Eq. (TTBTj) as 

n=0 V / 

= £ ( n I 2 ) 9i l+2) {{B ~ l) n+1 /i - (B + l) n+1 - B n+1 A 



n=0 \ 

= g[{B-l)fi-g[{B + 1)0-^3) A. (156) 

We consider the third term in the right-hand side of Eq. (I15ip in the following manner. 
At first, referring to Eq. (1112)) . we calculate the part including {R n+ i} in the third term 
of the right-hand side of Eq. (1151)) as 

( n + 3 \ >+3) p 
SI 2 J* R " +1 

= -£("t 2 )( n 2 4 )«!" +4, [i,B" +1 ]C 

oo 



+ei n V "t 5 Ur 5, u,B" +i a 



n=0 
oo 



-El n V T 9 r 6, [i,s" +i ]c 



n=0 

oo 



+ei "t 5 rt 7 Ur 7, [i,B" +i ]-4 



n=0 



(157) 



Here, we prepare the following formula: 



n + m + 1 \ I n + m + 3 \ ( n +™+3) n +i 



I m 2 M 



n=0 



1 1 d m ... . (m + l)(m + 2) (m+2 ) „ , s 

= 2^^W-- f ^ for- =1,2,3,.... (158) 

From Eq. ( 11021) and the above formula, we can rewrite Eq. ( 11571) as 
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1/1 d 2 



X 



2 V2! dx 2 
1/1 d 3 



2^31 dx 3 
1/1 d 4 



2 Ml dx 4 



g'({x 



A - 0) 

x =b-i dx 

. 1 d 2 



~ l~ g '^ x 
x =b+i ax 



x=B-l 



2\dx 2 
1 d 3 



x=B-l 



v 



x=B+l 



x=B-l 



fi 



3\dx 3 ' 
1 d 4 



x=B 

2\dx 2 
1 d 3 
3\dtf 
1 d 4 

x= b+i 4! dx 4 



A)C 



v 



x=B+l 



5( 

4( 



x=B-l 



d/dx It 



x=B-l 



4\dx 4 

g'[(B - 1) - e-^i(z) 
+ ^( J B-l) + e d / d V(^ 



x=B 



A) A 



x=B 



A)C 



x=B 



A) A 



x=B 



x=B 



9l 



KB)) 



-e~ d/dx g'{(x) 



x=B+l 



+ g'{{B + 1) + e~ d/dx g'{{x 



x=B 



-9i(B)m 



\(e d/dx g'l{x) 



x=B+l 



g'((B + 1) - e d/dx g'l(x 



x=B 



g'l{B))P 



l(g'{(B-2)-2g'((B-l)+g'{(B))(i 2 



+-{g'l(B - 1) - 2g'{(B) + g'l{B + 
+\{g'l{B-l)-2g'l{B)+g'l{B + l))vfi 
+ \{g'[{B)-2gl{B + l)+gl{B + 2)y. 



(159) 



Next, using Eqs. f)113D and (1122ft . we compute the part including {S n+ i} in the third 
term of the right-hand side of Eq. (11511) as 



E 

n=0 V 

oo 

-E 

n=0 



n + 3 \ ( n +3) 



g\ 



n + 3\ (n+3) 

o gi 



B - l) n+1 - (B + l) n+1 )D 



-Ug'{{B-l)-g'{{B+l))b. 



(160) 



Using Eqs. fl 1 2 j) and (11131) . we compute the fifth term in the right-hand side of 
Eq. (HSU as 



E 

n=0 



n + 4 \ ( n+4 ) 



n+l] 



P ( " 3 4 ) ~ 1)n+1/i " ^ + 1)W+lz> " 



x=B-l 



d 3 , , 

»-dx 39l{X ' 
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d 3 



x =B+i dx A 



x=B 



A). 



(161) 



Second, we examine the part which includes {[A, S n ]} in Eq. (I149p . From Eq. (11221) . 
after slightly tough calculations, we obtain 



[A,S n ] = 2 " \[A,B n - 1 }D + 2i 3 J [A, B n ~ 3 ]D + ... + 2 



n 

v 3 / 

+ ({B + l) n -(B-l) n )[A,D] 
for n = 0, 2, 4, ... (even), 



[A,£ n ] = 2 (J [i,B»- 1 ]£ + 2 (J \[A,B n - 3 ]D + ... + 2 



+ [(B + l) n -(B-l) n )[A,D] 
for n = 1, 3, 5, ... (odd). 

Here, we prepare the following formulae: 



n=0 

1 



(n+2m-l) [ n + 2m - 1 
I 2m - 1 



X 



2m- 1 



(2m- l)!cfe 2m - 



n 

n — 1 



n 

n - 2 



for m = 1, 2, 3, 



— + — — + — — + ... = -(e d/dx - e- d/dx ). 
dx 3! <ix 3 5! dx 5 2 



From Eqs. f |T02l) . f lT62D . ( 11631) and (HMD , we obtain 



[i,S 2 ]D 



(162) 



(163) 
(164) 



n=0 

oo 



^^((E+ir-^-irjtAD] 

n=0 

^n) ( n + l 

1 



71=0 



(n+3) / ri + 3 



71=0 



(n+5) / n + 5 



) [A,B n ]D 



n=0 



+ ... 



(^(B + l)-^-!))^] 



1 d 3 



x=B+l 



x=b-i dx 

1 d 3 

x =b-i 3! dx A 
1 

x =b-i 5! cfe 5 



0-4-9i(x) A)b 

dx X =B 1 

x =b+i 3! dx s 
1 

V ~ — — Ql{x) 



AW 



x=B 



x=B+l 



5! <ix 5 ' 



AID 



x=_B 



+ ... 
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9i(B + 1) — gi(B — 1) ) [A, D] 



[(e d/dx -e- d/dx ) gi (x) 



x=B-l 
x=B 1 



V 



x=B+l 



9l (B + l)- gi (B-l))[A,D] 
+ ((gx(B) - 9l (B - 2))fi - (gi(B + 2) - 9l (B))u 
-( gi {B + l)- gi {B-l))A)D. 



(165) 



Putting together Eqs. AIM]) . (|T54"1) . ( 1T55D . (|T56l) . ( TT59D . (TTHTB . ( 1T6TD and (lIBSjl . we 
obtain 



£^[i,i? n ] + [i,S n ]) 

n=0 



E 



-l) n / cT 



71=1 



n! ^c?x 



d n 

x=B~2 dx n 



+ "1 01 



(-l) n / tP 



n=l 
oo 

+ E 



. -2—g x (x) 

x=B-2 dx n 



d h 



n\ ^dx 

'^ri^dx^ 91 ^ 
+ ^(-l) n f d 



71=1 
°° 1 / cf ' 



f~; — tfiO^O 

\d,r n v y 



- 2- — 5fi(x 

a;=B-l as" 



x=B-l dx n 
d n 



x =b + dx n 



9i{x 
9i(x 



x=B 



fl 2 U 



x=B+l 



x=B-l 



d n , , 



71=1 

00 1 / cT 



n! ^c£r 
+ g (-1)™^ d n 



(-, — 9i( x ) 

\d.T n v y 



i=s-i dx' 



9i( x 
d n 



x=B ^ dx r 



x=B+l 

9i( x 



d n . . 
. -2-sH 
x =b-i dx n 



d n . N 

x =b dx n 



x=B+l 



x=B+l 

v\xi> 



VfJ, 



n\ ^dx 

E ~\{j~n9x{ 



71=1 

00 1 / d n 



(-, — 9i( x ) 



d n 

.~ 2 h 9 ^ 

X =B dx n 



d n , , 

x =b+i dx n 



n=l 
oo 



-2— gi (x 

x =b dx 11 



d n . s 

x =B+i dx n 



x=B+2 

v 3 



v 2 \x 



x=B+2 



n\ ^dx n 
+ ^6(l^ 9li 



71=1 

00 1 / d n 



71=1 

d 



cf' 

x=B-l dx 

d n 



:9i{x] 



-dx 

3\dx 3 
1 d 5 



x=B-l 
d 



dx n 



9i{x) 



x=B+l 

m 



Dp, 



5! dx 5 



9i(x 
9i(x 



x =b-i dx 

„ 1 d 3 , 



x=B+l 

d 



x=B+l dx 



A 



x=B 



x =b-i 3\dx 3 ' 

x=bJ-^ 9i(X] 



V 



x=B+l 



V 



x=B+l 



1 d 5 
5! dx 5 



9i{Xj 
9i( x ) 



A 



x=B 



A 



x=B 
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+...)[A,C] 

+ (g 1 (B + l)-g 1 (B-l))[A,D] 
+[(gi(B) - gi(B - 2))/i - ( 9l (B + 2) - 9l (B))u 
-( gi (B + l)- gi (B-l))A}D 
= (e- d ' dx gi {B - 2) - 9l (B - 2) - 2e- d ' dx g 1 {B - 1) + 2^(5 - 1) 
+e- d ' dx gi {B) - gi (B))jl 3 

+ (-e d ' dx 9l (B -2)+ gi (B-2) + 2e d/dx 9l {B - 1) - 2 9l [B - 1) 
-e d ' dx gi {B) + gi {B))fL 2 u 

+ (e- d / dx gi {B - 1) - - 1) - 2e- d ' dx g 1 {B) + 2«? 1 (5) 
+e- d / dx 9l (B + l)- 9l (B + l))jlvjl 
+ (-e d / dx 9l {B - 1) + - 1) + 2e d ' dx 9l {B) - 2 9l (B) 
-e d / dx gi (B + l)+ gi (B + l))fiu 2 

+ (e- d ' dx gi (B - 1) - - 1) - 2e- d ' dx g 1 (B) + 2< 7l (5) 
+e- d ' d *g l {B + l)- 9l {B + l))uft 2 
+ {-e d / dx gi {B - 1) + ^(S - 1) + 2e d ' dx gi {B) - 2 9l (B) 
-e d / dx gi (B + l)+ gi (B + l))ufiO 

+ (e- d/dx gi (B) - 9l (B) - 2e- d ' dx 9l {B + 1) + 2^(5 + 1) 
+e- <, / <fa 0i(S + 2) - 9l (B + 2))z> 2 /i 
+ (_ e d/«**0 1 ( J B) + gi (B) + 2e d l dx gx {B + 1) - 2^(5 + 1) 
-e^< ?1 (5 + 2) + <? 1 (5 + 2))i> 3 

+ (_ e -d/«fa^(B _ i) + gx (B - 1) + e~ d / dx 9l (B + 1) - ^(S + 1))l>£ 
+ (e d / dx 9l {B - 1) - ^(S - 1) - e d /<V(5 + 1) + + I)) I)v 
-f-(e d/dx - e- d/dx ) gi (B - l)fx - -( e d/dx - e - d/dx ) gi (B + l)v 

- X -(e dldx - e- d/dB )0i(S)i) [A, C] 
+ (^ 1 (S + l)- ffl (S-l))[i, J D] 
+ [(9i(B) - 9l (B - 2))/} - (^(S + 2) - ^(B))i> 
+ l)-< 7l (fl-l))i]£ 
= (- gi (B) + 3 ffl (S - 1) - 3 9l (B - 2) + ^(S - 3))/i 3 
+ (-<7i(£ + 1) + 3^(S) - 3 9l (B - 1) + ^(5 - 2))/i 2 z> 
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+ (-gx(B + 1) + 3 9l (B) - 3g x {B - 1) + g x {B - 2))fivfi 
+ (-gi(B + 2) + 3 9l (B + 1) - 3 9l (B) + 9l (B - l))/iz> 2 
+ (-gi(B + 1) + 3 9l (B) - 3 9l (B - 1) + 9l (B - 2))ufi 2 
+ (- 9l (B + 2) + 30i (B + 1) - 3^i (S) + 0i(S - 1))£>A^> 
+ (-0i (S + 2) + 3 9l (B + 1) - 3^i(S) + 9l (B - l))z> 2 /i 
+ (-0i (5 + 3) + 30! (5 + 2) - 30!(S + 1) + 9l (B))v 3 
+ (-0i (B - 2) + 0a(S) + 0i(S - 1) - 9l (B + 1))^ 
+ (</i(S) - 0i (S + 2) - 0!(S - 1) + 9l (B + l))£z> 
+2( 9l (B)- 9l (B-2))flb 
-2{g 1 (B + 2)-g 1 (B))uD 
-( 9l (B + l)- gi (B-lj)DA 

-( 9l (B + l)- 9l (B-l))AD. (166) 

Here, to compute the right-hand side of Eq. (I166p . we arrange fi in the left side of 
the product of operators and v in the right side of the product of operators. For the 
arrangement of operators, we carry out the following calculations: 



B n fi 3 


= V 3 (B + 3) n , 


B n li 2 v 


= fi 2 {B + 2) n z>, 


B n jj,vjx 


= fi 2 {B + 2) n v + jl(B + l) n D, 


B n liv 2 


= ft(B + l) n 2 , 


B n 0fi 2 


= A 2 (5 + 2) n z> + 2/i(5 + l) n (D + l), 


B n 0fi0 


= fi(B+l) n D 2 + B n Di>, 


B n 2 fi 


= fi(B + l) n v 2 + 2B n (t) + 1)0, 


B n Dfi 


= fi(B + l) n (b + 2), 


B n jl 


= KB+IT, 


B n AD 


= jl(B + l) n D-B n (D + 2)v, 


B n DA 


= fi{B + l) n (D + 2) - B n D0 



forn = 1,2,3,..., (167) 

vt) = (D + 2)P. (168) 

With making use of Eqs. (TT2Tj) and ([USD , substitution of Eqs. flTBTj) and (ITSgjl into 
Eq. ffT66|) yields 

oo 

= -(«|(a|E^! n) ([A^] + [i,5 n ])|a)|a) 

n=0 

= -8a 6 (a|(5|[-0i(S + 3) + 30i(B + 2) - 3 9l (B + 1) + 1 (5)]|a)|fi) 
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-2a 2 (a\(a\[- 9l (B + 2) + 3 9l (B + 1) - 3 9l (B) + 9l (B - l))D\a)\a) 
-Aa 2 (a\(a\[- 9l (B + 2) + 3 9l (B + 1) - 3 9l (B) + 9l (B - l))(D + l)|a)|d) 
-a 2 ( a \(a\[- 9l (B - 1) + ^(S + 1) + ^(fl) - 9l (B + 2)](L> + 2)|a)|a) 
-a 2 (a\(a\[ 9l (B) - 9l {B + 2) - 9l {B - 1) + ^(5 + l)]Z>|a)|d> 
-2a 2 {a\(a\[ 9l (B + l)- 9l (B-l))D\a)\a) 
+2a 2 (a\(a\[ 9l (B + 2) - <?i (£)](£> + 2)|a)|a) 
+a 2 (a\ (a\[ 9l (B + 2) - ^(B)]^|a)|a) 
-a 2 (a\ (a\[ 9l (B + 1) - ^(B - 1)](D + 2)|a)|d) 
+a 2 H(fi|[< 7l (£ + 2) - 9l (B)](D + 2)|a)|a) 
-a 3 (a|<5||>i(S + 1) - 9l (B - l)]D\a)\&) 
= 4« 4 (2« 2 + 3)gf ) (t) 

-12a 2 (2a 4 + a 2 -2)Q {2 \t) 
+4a 2 (6a 4 -3a 2 - 10)Q^ 1} (t) 

-4« 2 (2« 4 -3a 2 -4)gf } (t). (169) 

Similarly, we obtain Pg^it) as 

P$(t) = 4a 4 (2a 2 + 3)gf(t) 

-12a 2 (2a 4 + a 2 -2)g^ 2) (t) 
+4a 2 (6a 4 -3a 2 -10)g^(t) 

-4a 2 (2« 4 - 3a 2 - 4)gf (t). (170) 

Here, we pay attention to the following facts: In Eqs. H169[) and (I170p . replacing 
and Q^if) for n G {0, 1,2,3} with a certain constant that is not equal to zero, we can 
rewrite both of Pgiit) and P^it) as polynomials of a. Because all of the terms cancel 
each other out in these polynomials, we can finally confirm that they are equal to zero. 
This fact can be found in any order corrections of the perturbative expansion, such as the 
first, second and third order corrections. 
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